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ABSTRACT 


Maintenance, a key element of Integrated Logistics Support, plays a very 
vital role throughout an equipment/system planned life-cycle. Maintenance 
costs contribute a major portion of the life cycle costs of an equipment or 
system. Past historical records have shown that the cost associated with 
system maintenance is immense and usually takes up a large portion of the 
annual operating expenditure. Besides the costs, sound maintenance efforts 
contribute to better operational availability and reliability of a system. 
Therefore, the objective is to attain the proper balance of operations between 
performance and effectiveness, and logistics support, which largely includes 
maintenance, spares requirements, and the available budget. Adequate 
maintenance is essential to ensure the effective and economical support of an 
equipment or system. Therefore there is a need to design optimal maintenance 
policies to maximize appropriate measures of system effectiveness. These can 
be either to minimize operational and maintenance costs, to improve overall 
system reliability or to maximize operational availability. 

In this thesis, various maintenance scenarios are examined and the 
corresponding optimal maintenance actions are planned to take place at 
intervals chosen so as to maximize an appropriate measures of effectiveness. 
Preventive maintenance policies are also planned so that the overall reliability 
of the system is always kept above a specified minimum reliability level, while 
either keeping the cost per unit time to a minimum or maximizing the 
operational availability, subject to cost constraints. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this 
research may not have been exercised for all cases of interest. While every 
effort has been made, within the time available, to ensure that the programs 
are free of computational and logic errors, they cannot be considered validated. 
Any application of these programs without additional verification is at the risk 


of the user. 





I. INTRODUCTION 


A. MAINTENANCE 
Maintenance constitutes a series of actions to be taken to restore or retain 
an equipment or system in an effective operational state. [Ref. 1]. 


Maintenance occurs at three levels: 


e Organizational Maintenance. 
e Intermediate Maintenance. 


* Depot Maintenance. 


The Organizational Maintenance activity is performed by the ship's staff, 
for frigates or equivalent, using the onboard spares, tools, test equipment, and 
the documentation furnished. For smaller ships and patrol craft, 
Organizational Maintenance is carried out with the assistance of the Squadron 
Support Teams. Intermediate Maintenance is performed at fleet workshops, 
while Depot Maintenance is carried out at the naval dockyard or by 
commercial shipyards and firms. 

Any equipment or system introduced into the service needs to be 
maintained adequately so that it is readily available to perform a mission 


successfully at an acceptable performance level. Besides that, maintenance 


helps to extend the useful life of the equipment/system and ensures the safety 
of personnel using it. 


Basically, maintenance can be broadly divided into two types: 


* Preventive Maintenance. 


* Corrective Maintenance. 


Figure 1.1 shows the various maintenance forms. 


MAINTENANCE 


PREVENTIVE MAINTENANCE CORRECTIVE MAINTENANCE 


simple pianned minimal 
preventive overhaul repairs 
maintenance 





Figure 1.1 The Forms of Maintenance 


Preventive maintenance is scheduled maintenance that is organized and 
carried out in accordance with the documentation and records at a 
predetermined time following a predetermined plan. Preventive maintenance 


is normally associated with maintenance that occurs when the system is still 


operating. The main aim of preventive maintenance actions is to improve the 
overall system reliability and to avoid sudden or unexpected failures which 


may be catastrophic. Preventive maintenance can be sub-divided into: 


e Simple Preventive Maintenance. 


* Complete Overhaul/Preventive Replacement. 


Simple preventive maintenance actions usually consist of inspections, 
adjustments, tuning, cleaning, lubrication, minor calibration, and replacement 
of worn out components and parts before they actually fail. This type of 
maintenance is usually performed at the Organizational Maintenance level and 
often 1t does not affect the downtime of the equipment/system.(The downtime 
for these maintenance actions is usually negligible). 

Complete overhauls are carried out to bring the state of the deteriorating 
system back to "as good as new condition", and to prevent major impending 
failures. This type of maintenance is usually carried out at the Depot 
Maintenance level during planned refit periods and it usually incurs some 
significant downtime. Preventive replacements are carried out when the 
overhaul is not economical or when the parts are non-repairable. 

Corrective maintenance is the unscheduled maintenance carried out to 
restore a failed system to its operating state. Corrective maintenance can be 


subdivided into: 


e Minimal Repair. 


e Major Overhaul/Failure Replacement. 


Minimal repairs are minor repairs or component replacements carried out 
on failed components or assemblies that restore the system to its operational 
state without significantly improving the overall condition of the system. This 
type of maintenance is usually carried out at the Organizational Maintenance 
level. 

Major overhauls are carried out when the system experiences a sudden 
major failure and the work has to be carried out at the Intermediate 
Maintenance level or at the Depot Maintenance level. Failure replacement 
usually occurs when the system is beyond economical repair or when the parts 


are non-repairable. 


B. MAINTENANCE POLICIES 

Our main goal is to design a maintenance policy that will maximize 
appropriate measures of system effectiveness; this can be either to minimize 
operational and maintenance costs, to improve overall system reliability or to 
maximize operational availability (uptime). The most relevant measure of 
effectiveness used is to maximize a suitable measure of operational availability 
subject to budget constraints. Ideally, we would like to carry out many 
preventive maintenance and inspection routines to ensure that the system is 


in an optimum operating condition, to avoid sudden and catastrophic failures 


and damage to the system, and also to prevent accidents which may be 
detrimental to the people working in the vicinity. On the other hand, excessive 
preventive maintenance actions may not only be unnecessarily costly, requiring 
many manhours but may actually prematurely age equipment. Therefore, we 
need a balance or trade-off between two extremes. Our objective, then, is to 
select an optimal maintenance policy for the particular equipment or system, 
and to decide when to carry out the associated maintenance routines. Some 


conceivable maintenance policies are the following: 


* No preventive maintenance and corrective maintenance is carried out. 
The system is replaced at a fixed age. When the system fails before this 
age it is replaced with a new one. This type of maintenance is normally 
associated with low-level subsystems: items such as components, sealed 
modules or other non-repairable parts. Some common examples are the 
magnetron in the radar transmitter unit and the belts found in motors. 


e No preventive maintenance is carried out. The system is replaced at a 
fixed age. For failures that occur in between the planned replacement 
age, minimal repairs are carried out to restore the system to an 
operational state. This type of maintenance is usually associated with a 
system consisting of several components such that when a component 
fails the system fails, and replacing the component restores the system 
to operation. Minimal repairs are often carried out in the field, i.e. 
onboard ships or at a forward air base. 


e The system is renewed at a fixed age either by replacement or overhaul 
after which the system is "as good as new". Failures in between the 
planned age replacement can be classified into type I and type II. Type 
I failures are simple failures which are remedied by minimal repairs 
using the support elements onboard ships, and Type II failures are major 
failures which require base support facilities and are rectified by part 
replacement or complete overhaul. 


¢ Both preventive and corrective maintenance actions are employed. 
Preventive maintenance is planned at some time interval to improve the 


reliability of the system. Corrective maintenance actions are carried out 
whenever failures occur. 
Analysis of system availability and maintenance costs under the above 


options will be carried out subsequently. 


C. LIFE DISTRIBUTIONS 

Any system or equipment put into service will operate satisfactorily for 
some period of time and then cease to operate satisfactorily (fail). This 
operating time is not known in advance and is usually treated as a random 
variable. 

Let T be the period of satisfactory operation. We say T is the time to 
failure or life length of the system. T is a random variable and has a 
distribution F called the life distribution. 


Often, life distributions fall into three main categories: 


° No wear - such as the Exponential Distribution. 
° Wear out - such as the Weibull Distribution with shape 
parameter > 1 , Gamma Distribution with shape 
parameter > 1. 
. Wearin - such as the Weibull or Gamma Distribution with 
shape parameter « 1. 
Distributions which exhibit wear out are often modeled as JFR (increasing 


failure rate) distributions. Distributions which exhibit wear-in behaviour are 


often modeled as DFR (decreasing failure rate) distribution. 


, Other classes of distributions which exhibit wear out or wear-in such as 
NBU (new bettter than used) or UBN (used bettter than new) can also be used 
and are discussed in detail in Barlow and Proscham. [Ref. 2]. 

Most of the systems which are newly installed onboard the ships initially 
have a decreasing failure rate, sometimes referred to as the infant mortality 
phase or running-in period. Then the failure rate becomes rather constant for 
Some time and finally increasing, exhibiting wear out. In reliability, such 
failure rate functions are said to have a "bathtub" shape. The running-in 


period is usually under contractual obligation. 


D. FAILURE RATE 

The failure rate or hazard rate is one of the most important statistical 
characteristics of any equipment/system frequently used in maintenance or 
replacement studies. 

The failure rate, h(t), is usually defined as the "instantaneous" conditional 
probability of failure at age t, given that it has survived to age t. When F has 
density f, h(t) is given by 

Е (Е) 
BG) ся (1.1) 
1 - Е (+) 
The failure rate h(t) can also be expressed in terms of the hazard 


function, H(t) 


ан (+) 
h(t) » ———— (1.2) 
dt 


The hazard function is related to the distribution function by the 
following relationship 


H(t) + - log (1 - F(t)) (1.3) 


E. RELIABILITY 

When deciding upon system maintenance policies, the frequency of 
maintenance actions becomes a significant parameter. The frequency of 
maintenance for a given system is highly dependent on the reliability of that 
system. In general, as the reliability of a system increases, the frequency of 
maintenance actions will decrease; conversely, the frequency of maintenance 
actions will increase as system reliability is degraded. [Ref. 1] 


The reliability function or survival function R(t) 1s given by 
R(t) = 1 - F(t) = Г RO dt (1.4) 


R(t) is the probability that a new system will perform its mission 
satisfactorily for at least a certain time t. If T is the time to failure of the 
system then 

R(t) = P(T > t) (1.5) 

Sometimes we are interested in the chance of survival of the system in 

the future given that it has survived up till now. This is called the conditional 


survival function. So if T is the time to failure of the system with the survival 


function [1 - F(t)], then the reliability of the system, R(t;a), at some age a 2 0 
is given by 
(1 - F(tta) ] 
о (1.6) 
[1 - Е(а) ] 

Е. AVAILABILITY 

Availability is a measure of system readiness and it is one of the most 
important measures of effectiveness usually employed in mission-oriented 
Situations especially in the military environment. Operational availability is 
the probability that a system or equipment, when used under stated conditions 
in an actual operational environment, will operate satisfactorily when called 
upon [Ref 1]. System availability is influenced both by the inherent failure- 
proneness of the system and by the time and resources (support elements) it 
takes to restore a failed system to service. [Ref. 3]. Times to failure or ’up 
times’ and to restoration or ’down times’ may vary considerably, and not 
necessarily independently, depending upon the mode of failure, the time 
required to diagnose the failure, availability of special tools, test equipment, 
and spare parts, and the proper documentation and the required personnel 
skills. The long-run availability or steady state is expressed as follows: 

E[U] 


A, ШО (1.7) 
E[U] + E[D] 


where 
E[U] is the expected uptime of the system 
E[D] is the expected downtime of the system 
For a system operating at sea the availability at equation (1.6) is usually 


expressed as 


MTBM 
a en, (1.8) 
MTBM + MDT 


where 

MTBM is the mean time between maintenance 

MDT is the mean downtime, which includes the mean active 
maintenance time (M), expected logistics delay time (LDT) and the expected 
administrative delay time (ADT). The mean active maintenance time includes 


the expected time for preventive and corrective maintenance. 


G. MAINTENANCE TIME DISTRIBUTION 

The time required to carry out simple preventive maintenance actions and 
overhauls can generally be modeled as normally distributed with mean p and 
standard deviation o. Most of these tasks are standard and are carried out in 
accordance with the planned maintenance schedules; which stipulates the 
procedures to follow, spares, material,tools and test equipment that are 
required to perform the maintenance actions. [Ref 4]. The tasks usually 


requires a fixed amount of time to accomplish with very little variation. 
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The time required for corrective maintenance actions can be divided into 


three basic categories: 


* Active repair time. 
* Logistics delay. 


* Administrative delay. 


1. Active Repair Time 
Active repair time depends on the environment, state of equipment 
(hot or cold), and skill level of the technician; and it can be sub-divided into the 
following categories: 
e Recognition or detection ( often the time until actual occurrence of a 
failure and its recognition is not known). 
e Localization or diagnosis. 
e Correction or repair. 


e Verification or check. 


2. Logistics Delay Time 
Logistics delay time constitutes downtime that is expended while 
waiting for the availability of a spare part, waiting for a special tool or test 
equipment to perform repair, waiting for transportation, and waiting to use a 


facility required for the repair. 
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3. Administrative Delay Time 
Administrative delay time constitutes downtime of adminstrative 


nature, such as personnel assignment priority and organizational constraint. 


4. Distribution 

The distributions most commonly used to describe the downtime for 
corrective maintenance actions are exponential and log-normal. The 
exponential distribution tends to fit the type of equipment that requires 
relatively short durations of repair and usually corresponds to the replacement 
of a failed unit. Occasionally, much longer times may be required for major 
repair or for spares. The lognormal distribution is useful for situations where 
there are few downtimes of short duration, a large number of downtimes 
closely grouped about some modal value and a few downtimes of long 
durations. 

If X, the downtime, is a random variable having the lognormal 
distribution given by the probability density function 


_(n x - py 
SR) = 1 е 5 Osx<o (1.9) 


ху2ла? 


then the logarithm of X is Normal with mean p and variance o^. In the thesis 





We shall study several maintenance policies applicable to systems onboard 
ships. In Chapter II we introduce optimal maintenance for three types of 


policies. These are considered in detail in Chapters III, IV and V. 


II. OPTIMAL MAINTENANCE POLICIES AND MATHEMATICAL 
MODELS 


In the commercial environment determination of the optimum 
maintenance policy and time is of great economic importance. However, in 
many military situations, failure of a system in an operational environment is 
not only going to be more costly, but dangerous and may jeopardise the success 
of a mission. If a system has an increasing failure rate, such as the failure rate 
of a Weibull distribution with shape parameter > 1, it may be wise to replace 
or overhaul the system before it has aged too greatly. [Ref. 5]. This is very true 
for systems onboard the ships especially when they are operating many 
hundreds of miles from their home base and hence the support elements are 
not completely and readily available. A failure at sea may be catastrophic, in 
terms of cost and operational requirements. Although we can not completely 
avoid failures, however we can reduce the chance or probability of such 
catastrophic failures. This can be done by studying the failure distributions 
and then employing appropriate maintenance actions to maximize the various 
measures of effectiveness. 

The appropriate maintenance actions could be either preventive 
maintenance, failure replacement, minimal repair and preventive replacement 


or complete overhaul. The objective of a maintenance policy is to find a 
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sequence of times for carrying out the various maintenance actions that 
maximizes the appropriate measures of effectiveness over the operational and 
maintenance cycle of the system. 

The most commonly used maintenance policy is the policy based on age, 
usually referred as the age replacement policy. Sometimes the maintenance 
policy is based on the running hours of the system or equipment. If a system 
consists of many identical components, then the maintenance of these 
components are done in a block or group and is called the block replacement 
policy, such as the replacement of the diodes in the exciter unit of an 
alternator where the accessibility is poor. 

A recent survey of Preventive Maintenance Models for stochastically 
deteriorating single-unit system [Ref. 6] highlighted the use of some 
optimization models for repair and replacement policy evaluation. Most of 
these models were based on minimizing the long-run expected costs per unit 
time of replacement and minimal repair as the measure of effectiveness. The 
basic minimal repair model developed by Barlow and Hunter [Ref. 7] has been 
generalised and modified by many authors to fit more realistic situations. 


Minimal repair models generally assume the following: [Ref. 6]. 


* The system's failure rate function is increasing. 


* Minimal repairs do not affect the failure rate of the system. 
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e The cost of a minimal repair C, is less than the cost of replacing the 
entire system C,. 


e System failures are immediately detected. 


The long-run expected cost per unit time using a replacement age t for 


the basic model is given by 


— (2.1) 


СО 
t 


where N(t) represents the expected number of failures (minimal repairs) 
during the period (0,t]. 

Using the basic minimal repair model as developed by Barlow and Hunter 
[Ref. 7], Tilquin and Cléroux [Ref. 8] investigated an optimal replacement 
policy for the case where an adjustment cost C,(ik), incurred at age ik, 
1 = 1,2,3,... and k > 0, is added to the basic costs C, and C. They showed that 


the long-run expected cost per unit time is given by 


Gu S СО) (2.2) 
t 


C(t) 


where 
v(t) 


С” 009) = У са 


із0 


and v(t) represents the number of adjustments in the period (0,t]. Tilquin and 
Cléroux [Ref. 8] showed that the global minimum for equation (2.2) exists in 


the interval [0,29] when the life distribution is IFR. 
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In the thesis we will examine the various maintenance scenarios here 
called the policies and then based on these we will formulate the appropriate 
mathematical models using the stochastic and reliability theory. This 
mathematical models will depend on the desired measures of effectiveness 
required. From these models we can obtain the times for carrying out the 


appropriate maintenance actions. The measure of effectiveness that will be 


considered are: 


e Minimizing the costs. 
e Maximizing availability. 


e Mission reliability. 


The maintenance actions will be different for various measures of 
effectiveness and it is up to the Decision Maker to select which one is suitable 
for his scenario. 


The following maintenance policies are of interest: 


e Policy I Age Replacement. 
e Policy II Minimal Repair with Age Replacement. 
e Policy III Minimal Repair, Failure Replacement / 


Overhaul, Preventive Replacement / 
Overhaul. 


These policies are discussed in the following chapters. 
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III. POLICY I (AGE REPLACEMENT) 


The system is replaced at the time of failure or at some fixed time t, 
whichever comes first. (Instead of replacement, we could also overhaul the 
system,which on completion of overhaul is assumed as "good as new’). 

This type of scenario is usually associated with the repair by replacement 
policy often adopted at sea and applies to modules and sub-assemblies of 
equipment which requires support elements not available at sea. In lieu of 
repairing these modules and sub-assemblies at sea, sufficient spares are 
carried onboard or prepositioned at the forward operating areas so as to 
accomplish the respective missions successfully. The optimal number of spare 
requirements are based on the measures of effectiveness desired. 

The age at which the operating system is replaced depends on the 


following factors: 


* Failure distribution. 
* Costs of failure and preventive replacement. 
* Downtime of failure and preventive replacement. 
* Measure of effectiveness: 
- Minimize costs. 


- Maximize availability. 
- Mission reliability. 
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When evaluating this policy the following assumptions are made: 


Planned replacements are less costly than failure replacements. 


The mean downtime for a planned replacement is less than that for a 
failure replacement. An unexpected failure may incur additional Logistics 
Delay Time (LDT) and Administrative Delay Time (ADT), especially ifthe 
failures occured whilst the ship is at sea. 


The system exhibits an increasing failure rate distribution i.e. h(t) 
increases as t increases. 


The cost and downtime associated with simple preventive maintenance 
actions and minimal repairs is negligible. 


Preventive maintenance actions and minimal repairs do not improve the 
reliability of the system. 


Let T be the time to failure of the system. T is a random variable with an 


IFR distribution function F(t). A cycle is completed everytime a replacement 


is made. It can be either a failure replacement or a preventive replacement. 


(The system probabilistically starts over again and each replacement 


constitutes a renewal). By using the Renewal Reward Process, Ross [Ref. 9]: 


Expected long run average cost, C(t,), is given by 


Expected cost incurred during a cycle 
CC. (e). Te 
i Expected length of cycle 


Е [С] 
Y ===> (3.1) 
E [I] 
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where 
C is the total maintenance cost incurred during a cycle 
L is the length of a cycle 


Let t, be the planned replacement age 


failure preventive failure 
replacement replacement replacement 





{ 
C = ( (3.2) 
( 


where C, is the cost of preventive replacement 


C, is the cost of failure replacement 


С; > С, 
{ + + Dt 
Dee а ” (3.3) 
{ R+ T ТИНЫ 


where К, 15 the time of preventive replacement 
R, is the time of failure replacement 
R, > R, 

Then, 

C[1-FG,)]+CAFG,)] 


———— — — (3.4) 
[R +t 1 -FERFE + | даг 


са) = 
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and the above equation can be simplified to the standard form as shown in 


Barlow [Ref. 7] and Jardine [Ref. 10]. See Appendix F for c etails. 


C,[1-F(,)]+C{F()] 


C(t РР Во. 3.5 
мед RU -FED REG) +M) Си 

where M(t,) is the mean life during a cycle and is given by 
M(t,) = Ju -F(0))dt (3.6) 


The optimal preventive replacement age can be found by finding the value 
of t, that minimizes the cost function in equation (3.4) and the optimal value 


of t, is that value that satisfies the following equation [Appendix B] 
C, 
h(t,) f. ?[1-F(t)dt-F(t) = == * д1 = P Ж Гр РА "Өй 
C, C, 


The other measure of effectiveness that the Decision Maker is often 
interested is the availability (the probability the system is up at any time t). 

Assuming that there are only two states, that is the system 1s either up 
or down, then from the Regenerative Process Ross [Ref 9], Availability,A(t), is 
the Expected amount of time the system is up during a cycle divided by the 
Expected time of a cycle and is given by 


Mean life during a cycle 
A(t) = я 
Expected length of a cycle 
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M(t) 
RjJ1-F(t)] *RF(t) -M(t,) 


The optimal preventive replacement age can be found by finding the value 


A(t,) = (3.8) 


of t, that maximizes the availability function in equation (3.8) and the optimal 


value of pe is that value that satisfies the following equation. [Appendix В] 


Mm 


һа) 


The other measure of effectiveness of interest is the mission reliability 





M(t) = [ К, + Ға, 


(3.9) 
(R,-R) 
which is defined to be the probability that the system will complete a certain 


mission of duration d when it is at age t. This is given as follows 


R(t * d) 
R(t,d) = 


R(t) 
1 - F(t + d) 


а (3.10) 
2 - ЕСЕ) 


In military applications we usually like to maximize the availability of the 
system subject to some budget constraint C: 
max A(t,) 
sau Belt)» <s С (3.11) 


t, > 0 


P 


Some of the continuous life distributions that are commonly used to model 


the increasing failure rate of the system are: 
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a. Weibull. 

b. Gamma with shape parameter > 1. 

c. Log-Normal (depends on parameter). 

If we do not know the probability life distribution then we have to resort 
to a nonparametric approach [Ref. 10]. 

It should be noted that if the failure rate of the system is constant 1.e. it 
exhibits the Exponential distribution then by virtue ofits memoryless property 
we do not carry out any preventive replacement no matter what because a new 
system is just as bad and good as an old one. 

For our case let us assume that from historical data we know that this 
system has a Weibull qu s with shape parameter a and scale parameter 


à. Then the probability distribution function F(t) is given by 


-aD 
F(t)=1-e а>1 ,А>0 ,t20 (3.12) 


and the failure rate h(t) is given by 


aa-l 
h(t)=aA t а>1 ,А>0 „t20 (3.13) 


and the hazard function H(t) is given by 


H(t) = QD a>1, à>0, t20 (3.14) 
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It should be noted that if the shape parameter alpha = 1.0, then the 
failure rate, h(t) = A, which is a constant and this corresponds to an 
Exponential distribution whose distribution function is given by 

Fit) = 1 - е (9 1>0,120 (3.15) 

Figures 3.1 and 3.2 shows the plot of the Weibull density function and the 
corresponding hazard function respectively. It can be observed that as alpha 
increases the failure rate increases more rapidly. 

It should be noted that for a Weibull distribution 


воа. 
ма) - ЈА Pe dt (3.16) 


has no simple elementary closed form solution, but can be expressed in terms 


of the incomplete Gamma function. 


A. NUMERICAL ILLUSTRATION 


A case example was taken with the parameters as follow: 


C, = $ 25000 

С, = $ 37500 

R, = 8 hours 

R, = 16 hours 
a = 3.0 


À - 1/1390 hours 
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Figure 3.1 Weibull Density Function f(t) for a > 1, 
and 1/lambda = 1390 hours 
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Hazard Function of the Weibull distribution H(t), a > 1 
and 1/lambda = 1390 hours 


Figure 3.2 


24 


The above example was extracted from [Ref.12]. This example was based 
on the overhaul of aircraft engines. This same principle can be applied to the 
overhaul of motors, compressors pumps and weapon systems onboard ships. 

An IMSL subroutine DQDAG was available on the IBM3033P main frame 
computer at the Naval Postgraduate School, [Ref. 13], which could integrate 
the function very accurately and expeditiously by using a globally adaptive 
scheme based on Gauss-Kronrod rules. The estimate of the absolute value of 
the error using this scheme was 10’. A brief description of this subroutine is 
shown in Appendix C. 

The cost function and the availability function as in equations (3.5) and 
(3.8) respectively were plotted against time. Figures 3.3 and 3.4 shows these 
plots respectively. It is observed that the cost is minimum at about t. = 1453 
hours and at a rate of about $28.95 per hour. At this critical value of t, (1453 
hours) the availability is about 0.9883. It is also observed that the availability 
is maximum at t, = 1126 hours giving an availability of 0.9888 . For this 
maximum value of availability we need a budget of at least $29.92 per hour. 
Based on this information it is up to the Decision Maker to choose the optimal 
time t. to carry out preventive replacement or overhaul. 

From the graphs in Figures 3.3 and 3.4 it is observed that the plots are 
fairly flat and this gives some flexibility to the Decision Maker. The 


replacement age need not occur exactly at the optimal age. 
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Figure 3.3 Plot of the Cost Rate Function, C(t,). 
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Figure 3.4 Plot of the Availability Function, A(t,). 
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The optimal values of t, for the minimum cost and the maximum 
availability as mentioned above also satifies the equations (3.7) and (3.9) 
respectively. The values are 1453.45 and 1126.38 hours respectively. 

The Decision Maker may also be interested in the mission reliability of 
the system / equipment. The graph in Figure 3.5 shows the probability that the 
system / equipment will sustain 24 hours of continuous operation successfully 
when the equipment is at a certain age, say t hours. Suppose the Decision 
Maker wants a reliability of not less than 0.95; then the optimum replacement 


age t, is 1371 hours. 
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Figure 3.5 Plot of the Reliability Function, R(t). 


From the plot in Figure 3.5 it is observed that the reliability decreases as 


the equipment ages. This shows that even if the equipment does not fail and 
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if it has to sustain an operation for a specified duration of time, an older 
equipment will be less reliable and thus may not fulfill an operational mission 
successfully. 

Table 3.1 shows the values of the cost rate, availability and the reliability 
for the replacement age around the optimal values. It is observed that the 
availability value is very close whether the replacement age is at 900 hours or 
at 1600 hours. This is mainly attributable to the small values of the downtimes 
for the preventive replacement, R, and failure replacement, R,. 

The three measures of effectiveness as discussed above are important 
criteria to determine the amount of spares required to be carried onboard or 
prepositioned at forward operating areas. When a ship is assigned to an 
operational area, say for 2 months, away from the homeport, the ship will have 
to rely on the forward base for replenishment. The replenishment cycle occurs 
every two or three weeks and usually takes about two days. So if the optimal 
replacement age is close (not necessarily within) to these stand-off periods we 
can undertake the replacement actions during these periods. It should be noted 
that if we do not carry the replacement action when it is due then if we 
continue to delay these actions everytime they are due, on the long run we are 
going to experience high costs, poor availability and low reliability. 

Since we have three measures of effectiveness and the respective optimal 
replacement ages, it is up to the Decision Maker to decide which one he is 


going to give top priority. 
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TABLE 3.1 VALUES OF COST/HR, AVAILABILITY AND 
RELIABILITY FOR VALUES OF REPLACEMENT AGE CLOSE TO 
THE OPTIMUM. 


Replacement 
Age (t,) 


Availability Reliability 
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B. SENSITIVITY ANALYSIS 

Since we do not actually know the failure parameters, it is also important 
to do sensitivity analysis and see how the various measures of effectiveness 
vary with changes in those parameters. This will enable us to determine the 
parameters that we need to estimate very accurately. 

We have assumed that the failure distribution is Weibull. Therefore, the 
parameters of interest are the shape parameter a and the scale parameter A. 
Besides these, other parameters are the expected costs and the downtime of 
preventive and failure replacement. In order to see the changes in the 
measures of effectiveness we hold all the other variables constant and only 
vary the parameter of interest. It would be possible to use Experimental 
Design techniques such as the Factorial Designs to study the main effects on 
the respective measures of effectiveness. This may help us to identify the 


parameters that are more sensitive and shall be estimated accurately. 


1. Alpha 
Table 3.2 shows the optimal value of t, for the minimum cost rate, 
the maximum availability and the reliability of 0.95 for sustaining a mission 
of duration 24 hours at age t; and the respective optimal measure of 
effectiveness for various values of the shape parameter alpha. The graphical 
plots are shown in Figures 3.6, 3.7 and 3.8 for the cost rate, availability and 


the reliability functions respectively. It is observed that as the value of a is 
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TABLE 3.2 OPTIMAL REPLACEMENT AGE FOR VARIOUS 
VALUES OF ALPHA 


m 


Шо 15260 | 29.23 | 11566 | 0.9886 | 14245 


a 14535 | 2895 | 11264 | 09882 | 13711 a 


3.5 1340.7 28.24 1083.4 0.9893 1289.7 





increased the replacement age gets shorter in order to minimize the costs and 
maximize the availability. This is also obvious from Figure 3.2 where the 
hazard increases as alpha increases and as such the failure rate also increases, 
which in turns requires replacement action early so as to optimize the 
respective measures of effectiveness. From Figure 3.7 it is observed that the 
curves are much flatter for low values of a and for a = 3.5 the curve falls quite 
rapidly on both sides of the optimal value. Figure 3.8 shows that if we want a 
minimum reliability of 0.974 then the optimal replacement age is 
approximately 980 hours and at this value of reliability it is insensitive to the 


shape parameter alpha. 


31 


Cost ($/hr) 


1 
з ‘ L ' 
‘ 
> > > ES 2 S > 
cr RAI rr ar rr rro ШЧ о МИ: ро сос зе сорс ес то 
1 1 * , 1 


i 2e 


a 
+ 
2 
g 
g 
> 
< 


800 1000 1200 1400 4600 1800 2000 
Time Chr) 


Figure 3.7 Plot of the Availability Function for Values of o. 
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Figure 3.8 Plot of the Reliability Function for Values of Alpha. 


From Table 3.2, it is observed that if the true value of the shape 
parameter & was 3.2 and in the analysis we estimated it to be 3.0, then we 
would have lost an availability of 0.0002 and would have incurred additional 
cost of $0.28 per hour. On the other, if we have estimated it to be 3.5 then we 
would have gained an availability of 0.0003 at a lower cost by $0.43 per hour. 
Therefore it is much better to estimate the shape parameter higher than lower 


and hence carry out the replacement or overhaul action earlier than later. 


2. Lambda 
Table 3.3 shows the optimal value of t, and the respective optimal 


values for the cost, availability and 0.95 reliability to sustain a mission of 24 
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hours duration. The graphical plots are shown in Figures 3.9, 3.10 and 3.11 for 
the cost rate, availability and the reliability functions respectively. 


TABLE 3.3 OPTIMAL REPLACEMENT AGE FOR VARIOUS VALUES 
OF LAMBDA (NOTE: MU = 1/LAMBDA) 


AVAILABILITY RELIABILITY 


1390 1422.2 29.58 1102.1 0.9885 1326.7 e 
1380 29.16 1118.3 0.9887 1356.3 


1390 | 1453.5 | 28.95 | 1126.4 | 0.9885 | 13711 
1463.9 | 28.75 | 1134.5 | 0.9889 | 1386.1 


1420 1484.7 | | 1150.7 0.9890 1416.2 





From th. above table it is observed that the optimal replacement age 
for achieving a reliability of 0.95 when subjected to a mission of duration of 24 
hours increases as p increases. This is because the higher the value of p, the 
reliability becomes better. Also when the value of n increases the optimal 
replacement age becomes longer for both the cost rate and the availability. The 
higher the value of p the failure rate for the Weibull distribution decreases and 
therfore the optimal replacement age becomes longer. A small difference in the 


value of p do not affect the replacement age drastically. However this is not 
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Figure 3.9 Plot of the Cost Rate Function for Values of n. 
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Figure 3.10 Plot of the Availability Function for Values of p. 
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Figure 3.11 Plot of the Reliability Function for Values of p. 





true for the case of alpha. Therefore, the shape parameter a need to be 


estimated more accurately than the scale parameter 7. 


3. Costs 
The costs for the failure replacement and the preventive replacement 
only affects the cost rate function. These costs estimates are quite easily 
available and they are usually fairly accurate based on previous data. Since we 
have not discounted costs and cash flow problems in the future it is necessary 
to observe the effects of the costs on the cost rate function. It is assumed that 
the increase in costs are proportional such that the ratio of replacement costs, 


C/C,, is the same. It is observed that the replacement age does not vary as 
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long as the ratio of the replacement costs, C/C,, remains the same. Equation 
(3.7) clearly shows that if the ratio of C/C,, is the same, there is no effect on 
the optimal replacement age bs However the cost per unit time increases 


proportionally as the replacement costs increases. 


4. Downtime 

The availability function greatly depends on the expected downtime 
of the failure replacement (R,) and preventive replacement (R,) assuming the 
parameters for the Weibull distribution are estimated accurately. In our 
example we have taken R, = 16 hours and the R, = 8 hours. Because these 
figures are small compared to the optimal replacement age the availability is 
very high. In reality a failure at sea may take days to repair taking all the 
logistics and administrative requirements into considerations, especially when 
the ship is not accompanied by any auxilliary vessel or support ship. That is 
the reason all vital equipments onboard have redundancies incorporated. We 
can carry spares onboard but may not have other support elements to rectify 
the defect. Usually the downtime for a preventive replacement is constant 
because it is done at the base with all the support elements and the work is 
repetitive following standard procedures. Therefore, in our analysis we shall 
only vary R, to see its effect on the availability. Table 3.4 shows the values of 
the availability if we carry out the preventive replacement at 1126 hours. It 
is observed that the availability of the system drops as R, increases, but is still 


quite high. 
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TABLE 3.4 AVAILABILITY FOR VARIOUS VALUES OF 
R, AT T, - 1126. 


Expected Downtime — 





C. EFFECTS OF SIMPLE PREVENTIVE MAINTENANCE 

In the above policy it was assumed that simple preventive maintenance 
actions did not improve the operational reliability of the system. In real life it 
is evident that these preventive maintenance actions may definitely enhance 
the system condition and, hence the reliability, but it will not restore the 
system to the original state i.e. to be as good as new. This was also emphasised 
in [Ref. 14]. However in reality there are situations for which preventive 
maintenance actions can degrade the system, that is by imperfect repair. Here 
we will concentrate on those actions that will improve the condition by a 
certain factor. So now we shall incorporate these preventive maintenance 
actions into our model. The questions to be asked is that when should we do 


these type of maintenance? How do we determine the improvement factor? 
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We can schedule the simple preventive maintenance actions whenever the 
operational reliability hits a minimum acceptable limit R,,, say 0.75 and then 
as mentioned above our reliability will improve by a certain factor, thus 
making the system ’younger’ but not necessarily new. Or we could also carry 
out a planned replacement when the operational reliability hits the Rmin value 
and bring the condition of the system back to new. The appropriate 
maintenance action to be taken at these critical times depends on the measure 
of effectiveness required. 

[Ref. 15] solved this problem by minimizing the expected cost rate as the 
measure of effectiveness. The number of simple preventive maintenance 
actions before a planned replacement action which gives a minimum cost rate 
is then obtained by evaluating the cost rate whenever the system reaches Rmin- 
Here we shall carry out the appropriate maintenance actions so as to maximize 
the availability of the system. 


Notation is as follows: 


1P simple preventive maintenance 
2P planned preventive replacement 
2C failure replacement 

T time to failure 


R(t.) probability of no failure during (t, „t,) 
Ft) probability of a failure during (t, ,,t.) 


F(t) probability distribution function of the time to failure 
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f(t) probability density function of the time to failure 


h(t) failure rate function of the time to failure 


К, time required for each maintenance action i=1,2 j=P,C 
С» cost of each maintenance action 1=1,2 ;j=P,C 

Ber minimum acceptable operational reliability limit 

t, ip time the system reaches Rmin 


Let I be the improvement factor, then if maintenance type 1P is carried 


out at time t, then the system age is reduced from t,, to ( t,-I ) and therefore 


t: = t, + t, (1 - I) 
and t, = t, * (t, - t) (1 - I) 
in general t, = t, + (1 — I)*"”` t, (3.18) 


As I — 1, the effect on the system approaches "bad as old" and 
As 1 > 0, the effect approaches "good as new" [Ref. 15]. 


The value of t, can be computed as follows: 
1 
L - = І (3.19) 
and from equation (3.19) we can compute t, fori = 2,8.,...,n. 
The value of I can be estimated from past records from the data collected 
on similar equipment where performance measurements are taken before and 
after the simple preventive maintenance actions. In some cases performance 


measurement techniques such as condition monitoring by vibration analysis 


are employed. This is largely used in rotating machinery. As mentioned earlier, 
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simple preventive maintenance actions constitutes cleaning, adjustments, 
replacing worn out parts, tuning and minor calibration. All these actions 
normally improve the reliability which is a function of the failure rate. So 
these actions may prolong the life of the system, and enable the system to be 
‘younger’. From experience, it is found that as the system becomes older these 
simple preventive maintenance actions becomes ineffective and as such the 
system may require a complete overhaul which on completion is as good as 
new. For our case here we shall consider that on completion of each simple 
preventive maintenance action, the system age is reduced from t, to t; - (t; - t;.,) 
(1 - I) where t; is the time of the simple preventive maintenance action. 

A cycle is completed each time a replacement takes place which could be 
either due to a failure replacement or a planned replacement. On completion 


of the replacement action the cycle probabilistically starts all over. 


1. Casel 


T >t, (Failure occurs after the type 2P maintenance) 


LP 1P 2P 
— 2 
0 t, t, ae t. 


The costs associated are: 
a) Expected costs of type 1P maintenance. 
If planned replacement takes place at t,, then there will be (n - 1) 


maintenance actions of type 1P and each costs C,p. 
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Then total cost = (n - 1) Ср 

But this event will occur if there are no failures in the first (n-1) 
intervals. The Probability of this event is = R,(t,)R,(t,)...R,(t,). It is assumed 
that the failures are statistically independent. The failure in each interval is 
independent from the failure in next interval. 


Therefore expected costs of type 1P maintenance = 
(n - 1): Ce [| RE (3.20) 
i=1 


where R;(t,) = 1 - Probability of a failure in the (n-1) intervals and is given by 
R(t) 1 - [F(t) - F(t; - (t; - t, )) (3.21) 

b) Expected costs of type 2P maintenance. 
Since there will be only one Type 2P maintenance in an interval, the 


expected costs of type 2P maintenance - 
Ср“ Пк) (3.22) 
i=] 


The expected downtime can be obtained by substituting the cost by 
downtime in equations (3.19) and (3.20). 
The downtime associated are: 


a) Expected downtime of type 1P maintenance = 


(n - 1) В," А) (3.23) 


i=1 


b) Expected downtime oí type 2P maintenance = 
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n 


Rp 1] қа) 324) 


і=1 


2. Case 2 
хат < t, 
1Р 1Р 1Р 2C 
0 t, t, m tua t, 


The costs associated are: 
a. Expected cost of type 1P maintenance 
If the failure occurs between t,, and t,, then there will be (n-1) 
preventive maintenance actions and the Probability of this event is 
= Probability of no failures in the first (n-1) intervals multiplied by the 
Probability of a failure between t, , and t,. Again we assume that failures are 
independent. Then the expected cost of type 1P maintenance - 
n i-1 
> [i - 1) ° Cip П R (t F (t)] (3.25) 


b. Expected costs of type 2C maintenance = 
п 1-1 
УГ С» TII 60) · Е) (3.26) 
i=1 j=l 


The downtime associated are: 


a. Expected downtime of type 1P maintenance = 
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n i-1 
$9 I6 - 9-5“ ЈЈ АС) : Еко (3.27) 
i=1 j=1 

b. Expected downtime of type 2C maintenance = 


n i-1 
УМА TI RE) Fe) (3.28) 
i=] j=1 


The expected operation time until t, is the sum of the average 
operation time until t„, and the average operation time during (t,,,,t,). 

The expected downtime is the sum of all the downtimes in equations 
(3.23), (3.24), (3.27) and (3.28). 

From our earlier discussion we know that the mean life of the 


system, assuming only 1 interval, is given by 


Mtt) - [t1 - Fo) d (3.29) 


and the above equation can also be written as 


Mt) = [rf dt + Вера. (3.30) 


and for n intervals the equation can be written as 


i-1 
>| к) if, NERO de Rap, - a (3.31) 
ізі 17 9475-1 


i-1 


Total Operation Time 
Now Availability = 


Total Operation Time + Total Downtime 


It will be shown later that if the expected downtime for a type 1P 
maintenance (simple preventive maintenance) is small then the availability 
function given by the above equation is a good approximation. 

The expected costs is the sum of all the costs in equations (3.20), 
(3.22), (3.25) and (3.26). Then the cost rate function, C(t,), is given by 


Total Expected Costs 
ES таннин — —— 
Total Operation Time + Total Downtime 


3. Numerical Illustration 
Now let us take a case example with the following data: 


Let the time to failure, T, follow a Weibull distribution with shape 


parameter a and scale parameter À. 


Let 
a = 3.0 
À = 1/1390 
Rp = 1 hour 
E = 8 hours 
R,c = 48 hours 


Using the above data and equations, a Fortran program was written 
[Annex A] to determine the optimum times to carry out simple preventive 


maintenance actions and complete overhauls or preventive replacement so as 
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to maximize the availability. The results are tabulated in Tables 3.5, 3.6, and 
3.7 for values of I = 0.1, 0.5 and 0.9 respectively. 
When I = 0.1, carry out simple preventive maintenance actions at 


times: 


e 917.6 hours 
e 1743.4 hours 


e 2486.7 hours 


and preventive replacement or complete overhaul at time 3155.6 hours to 
achieve a maximum availability of 0.9836. 

When I = 0.5, carry out simple preventive maintenance action at time 
917.6 hours and a preventive replacement at time 1376.4 hours to achieve a 
maximum availability of 0.9802. 

When I = 0.9. do not carry out any simple preventive maintenance 
action and at time 917.6 hours replace or carry a complete overhaul of the 
system so as to achieve a maximum availability of 0.9794. 

It is observed that as the improvement factor, I — 0, the system 
becomes "good as new’ after each simple preventive maintenance action and 
as the improvement factor,I > 1, the system becomes "bad as old” after each 


preventive maintenance action. 
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TABLE 3.5 REPLACEMENT TIMES FOR IMPROVEMENT FACTOR, 
I = 0.1 AND R;, = 8 HOURS, R,, = 48 HOURS, R,p = 1 HOUR, a = 3.0, 
p = 1390 HOURS. 


A] o 
а [шш | owe 
NENNEN | 
| и | я 
Cos oom | ome 
ШЕН БЕСТЕН ШЕГІН 


0.983206 


5621.0 0.983075 
5976.5 0.982955 


6296.5 0.982848 
6584.4 0.982751 
6843.6 0.982664 


14 7076.8 0.982587 | 


15 7286.7 0.982518 











10 
11 
12 

3 


1 
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TABLE 3.6 REPLACEMENT TIMES FOR IMPROVEMENT FACTOR, 
I z 0.5 AND R;, - 8 HOURS, R4, = 48 HOURS, R,p = 1 HOUR, a = 3.0, 
п = 1390 HOURS. 


Is 
ГГ = [= 
ШЕН БЕСТЕН ШЕГІН 
ШЕКЕ БЕСТЕН чен — 
БЕН БЕСТЕН БЕСІН 
2 | ән | чен ` 
кишш ишик 
а вне мно 


1835.1 0.978964 


1835.1 0.978938 
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TABLE 3.7 REPLACEMENT TIMES FOR IMPROVEMENT FACTOR, 
I = 0.9 AND R,r = 8 HOURS, R,. = 48 HOURS, R,, - 1 HOUR, o = 3.0, 
p = 1390 HOURS. 


us | avanan 
а | ma | oras | 
[s |77 | eme 
ШЕН ШЕКТЕН БЕКІ; ЖЕН 
Се | m ење | 


T" 1019.6 0.976426 
реж 1019.6 0.976123 
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IV. POLICY II (MINIMAL REPAIR WITH AGE REPLACEMENT) 


As an alternative assume that the system consists of several components 
and the system fails when one of the components fail, that is components are 
connected in series. We shall assume the components have independent IFR 
distribution. Then the system lifetime is also IFR. (The IFR property is 
maintained under the formation of the series system). Preventive replacement 
is carried out when the system reaches the age t,. Between the preventive 
maintenance, failures are repaired as quickly as possible (minimal repairs) 
either by replacing the failed component with a new one or repairing the failed 
part. We also assume that the system failure rate is not disturbed on 
completion of the minimal repair i.e. the failure rate is the same as before 
failure. If the failure occurs at time t < t, the failure rate of the system just 
after the minimal repair is h(t). Failures are detected immediately. 

This idea of minimal repair was first introduced by Barlow and it is 
described in [Ref. 7]. After that, many authors expanded on these ideas and 
formulated various models associated with minimal repairs; they used the 
expected cost rate as the measure of effectiveness. In this policy we will 
emphasize the availability of the system as our measure of effectiveness. 


However we will also formulate the expected cost rate models. 
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Let 0 x t be the system operating time since last replacement. Then the 
probability of a failure occurring in [t,t--dt] is h(t)dt where h(t) is the failure 
rate of the system. 

Let N(t,) = number of failures occuring in time (0,t,) 

When repairs are minimal, { N(t,) , t, 2 0 } is a Nonhomogeneous Poisson 
Process with intensity function h(t,). 

The probability that the system will experience n failures in the interval 
(0,t,) is given by 


Impr en? 


(4.1) 
n! 


P{ Ма) =n} 
where 

m(t) = EL NG) Ј? коа (4.2) 
m(t,) is called the mean value function. 
A. COST MODEL 


The age t, at which the operating system is replaced depends on the 


following factors: 


* Failure distribution (IFR). 
* Costs of minimal repair and preventive replacement. 


* Downtime of minimal repair and preventive replacement. 


o1 


° Measure of effectiveness: 
- Minimize costs. 
- Maximize availability. 
We have assumed that there are no simple preventive maintenance 
actions. 
-8,- све -R -R =—-R——- 


————a, T q IA 
t, 


x----failure 


Let C, be the mean cost of each minimal repair. 
C, be the mean cost of preventive replacement. 
R, be the mean time of preventive replacement. 
R, be the mean time of a minimal repair. 
and t, be the planned replacement age. 
C,>C, and R>R, 
Then, assuming that the cost of downtime is negligible, the total cost is 
given by 
ма) 
С- >, C, + С, (4.3) 
where C, is the cost of the i-th minimal repair and the length of the cycle is 
given by 


L = t, + R, (4.4) 
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By using equation (3.1) and assuming that the time for a minimal repair 
is very small compared to that of the length of a cycle, the expected cost rate 


function is given by 


c, [* ^ à «c, 
t, + É, 


and we wish to minimize the cost per unit time, so we set the derivative of 


Са) = (4.5) 


equation (4.5) to 0 and we obtain 
G 
вари, + Ва - f? hO d = E (4.6) 
f 


Then the value of te that satisfies the above equation is the optimal 
replacement age. 

For a time to failure, T, following a Weibull distribution with shape 
parameter о апа scale parameter X the expected cost rate function is given as 


below 


а 
са) - чю 10 (4.7) 
р* № 
For given values of C, , C. , R. , R, and the parameters a and A we can 
find the optimal values of t, that minimize the expected cost rate. 
Now let us take a case example with the data as follows, but first without 
preventive maintenance action: 


C, = $ 25000 
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R, = 8 hours 
R, = 1 hour 
a = 3.0 


À = 1/1390 hours 
Figure 4.1 is a plot of the cost rate function as in equation (4.7) and the 
optimal value of t, that minimizes the cost rate function is 3222 hours at a 


cost of $11.60/hour. 
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Figure 4.1 Plot of the Cost Rate Function. 


Table 4.1 shows the values of the cost rate for the replacement age t, 
close to the optimal value. From Figure 4.1 it is observed that the curve is 
fairly flat near the optimal value thus giving some flexbility to the decision 


maker. 


54 


oe “п — 


TABLE 4.1 VALUES OF COST/HOUR FOR VARIOUS 
REPLACEMENT AGE. 


Replacement Age 


(hours) 





Now suppose there are s identical components which operate and fail 
independently and have the same failure distributions with the same 
parameters. All of these components are replaced at time t, at a cost of C, . 


When each of the component fails it is replaced individually or undergoes 
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minimal repair without affecting the system failure rate as a whole. Then the 


expected cost rate function becomes 


5 С, Ji h(t) dr + C, 
t, + R, 


Now let us imagine that we carry out a simple preventive maintenance 


Cee (4.8) 


action at some time t, and after N of these preventive maintenance actions are 
carried out we either replace the system or carry out a complete overhaul 
which on completion is like new. Any failures between the preventive 
replacement are treated as minimal failures and are repaired quickly in 
negligible downtime at a cost of C, The cost of each simple preventive 
maintenance action is C,,. It is assumed that after each simple preventive 
maintenance action the system improves by a certain factor, so for simplicity 
we say that the system becomes ’younger’. Nakagawa [Ref. 16]. This means 
that if t is the time of a preventive maintenance action the failure rate on 
completion is h(t-x) where x is the amount of time by which the system has 
become ’younger’. However the failure rate after a minimal repair stays the 
same. The value of x can perhaps be determined from past historical records 
using performance measurements techniques such as condition monitoring by 
vibration analysis or some output parameter. However a methodology for 
characterizing the effective age reduction, x, remains to be developed. 
According to the above assumption the replacement age is Nt,, and 


there are N-1 simple preventive maintenance actions. 
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A cycle is completed at (Nt, + R,) and the process probabilistically 


starts all over. 


ee O 
0 ti t; Es Nt. 

f ---- minimal failure 

R, ---- expected time of preventive replacement 

t, 1 = 1,2,...,N-1 ---- simple preventive maintenance 

Nt ---- preventive replacement 


P 


( Note: At t, the age of the system is (t, - x) ) 
The total expected costs, E[C], incurred in a cycle is given by 


t, *(t, -x) 
x 


ЕС] - с, |?қфй > с.» с, | мда» с...» с, (49) 


for O < x <t, 


and this equation can be simplified to 
да t +(t -xJ 
ЕС] = С, > Іі a h(t) de + (N-DC, +C, (4.10) 


for 0 <x <t. 
and then the expected cost rate function, C(t,), is given by 
E[C] 
СЕ) 
ЕТІ.) 


where E[L] = Nt, + R, 
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For a time to failure, T, following a Weibull distribution with shape 
parameter a and scale parameter à ne expected cost rate function is given as 


follows 


N-1 
t, *(t, -x)i a.a-1 3 
CY [t av dk + (N-DC,, «C, 


са) - AES e E (4.11) 
P 


When C,, = 0 and x = 0, that is there is no preventive maintenance 
action, then the above equation is the same as equation (4.7). 

Now we shall plan to carry out a preventive maintenance at some 
time t, at a cost of C,, = $ 500.00. On completion of this simple preventive 
maintenance action, the age of the system becomes younger by some value x. 
When x = 0, it is observed that there is no improvement to the optimal planned 
replacement age but we have incurred an additional cost at $11.75/hour 
instead of the original cost at $11.60/hour. As x increases from 0, the optimal 
planned replacement age increases and also the cost rate reduces. This is for 
the case of carrying out only one preventive maintenance in between the 
planned replacement or overhaul age. We could carry on the same analysis 
incorporating more preventive maintenance actions and the results will give 
more improvement. However this greatly depends on the parameters such as 
the value of x and the cost of preventive maintenance assuming the other 


parameters remain fixed. 


o8 


In the real scenario as the system gets older there is more 
deterioration and on completion of each preventive maintenance action the 
value of x probably reduces over time; that is; it is some function of t, and 
there comes a stage at which preventive maintenance action will not improve 
the system any longer. For simplicity, however we assume that x is constant. 
It should be noted that there will be values of age for which (age - x) « 0 , we 
therefore take the maximum (0, (age - x)) as the age of the system on 
completion of a simple preventive maintenance action. 

The Figure 4.2 and Table 4.2 illustrate that simple preventive 
maintenance actions increases the optimal replacement age. 

Table 4.2 shows the time for a preventive replacement (TPR) and the 


time for a simple preventive maintenance (TPM) for various values of x. 


B. AVAILABILITY MODEL 
Now let us take the availability as our measure of effectiveness. From 
earlier discussion Availability is given by 
Mean Life During A Cycle 


A(t,) = 
Total length of cycle 


г, - , Ј" h(s) ds 


(4.12) 
t, г R, 


Ай) - 


апа we wish to maximize the availability function, so we set the derivative of 


equation (4.12) to 0 and we obtain 


59 


tax 
i 
Э 
© 
fs 
мч 
A 
ы 
к 
8 


3000 3500 4000 
Time Chour) 





Figure 4.2 Plot of the Cost Function with Preventive Replacement 
for Various Values of x. C,, = $500.00. 


TABLE 4.2 OPTIMAL PREVENTIVE REPLACEMENT AGE WITH 


PREVENTIVE MAINTENANCE FOR VARIOUS VALUES OF X. 
C,» - $500.00. 


"COSE 
(hours) (hours x 1) (hours x 0.5) 2 


—— зз И 
[om o | w [ox 
mr 
| 8 







3573 1787 
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h(t) [t, * RJ - |, h(t) dt = z (4.13) 


Then the value of t, that satisfies the above equation is the optimal 
replacement age. 
For a time to failure, T, following a Weibull distribution with shape 


parameter a and scale parameter A the availability function is given as follows 


а 
A(t) = t, ; к=; (4.14) 
р р 

The availability function in equation (4.12) is an approximation because 
we do not take into account of the downtime during a minimal repair, which 
in reality is not entirely valid. In the equation we assume that the system can 
still fail when the system is down (we are integrating over (0,t,) and some of 
that is downtime). However, this approximation is quite accurate if the 


downtimes are small, as is likely to be true in practice. 


We can find an exact solution to equation (4.12) as follows: 


Let 
A(t) = Probability the system is available at time t (age t) 
following the last replacement. 
h(t) = failure rate at time t. 
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p(t) = repair rate at time t. For our case this is assumed to be 
constant and denoted by pn. Then n = R,’. 
(It should be noted that it may take longer to repair an older system) 
Then, 
A(t+dt) = A(t)(1 - h(t)dt) + (1 - A(t))u(t)dt + o(dt) (4.15) 
where o(dt) represents higher order terms which are = 0 
From equation (4.15), arranging the terms we get 


Alt+dt) - A(t) 
————— = -h(t) A(t) + w(t) {1 - A(t)} (4.16) 








dt 
dA (t) 
= -h(t) A(t) - u(t) A(t) + p(t) (4.17) 
dt 
dA (t) 
+ {h(t) + w(t)} A(t) = u(t) (4.18) 
dt 
PELO Bl 09 * 1014 (4.19) 
fs the) + ap, t fo. (t) + wep) a, (4.20) 
:. A(De -А(0) = jr u(t,)e dt, i 
-f ihe) + -|авер + вар 
PAY ја, [ъс Така) + вера de (4.21) 


For a time to failure, T, following a Weibull Distribution with shape parameter 


a and scale parameter A, the failure rate h(t) is given by 


h(t) » aA*t* a>1,A>0,t20 (4.22) 
then subtituting A(0) = 1 and , p(t) = (R2! ; A(t) can be simplified as follows: 
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t b 


КА“ + = Ад“ - (Anyi +t - 2) 
Ай) =e A e КК, (423) 
R, 0 
-[an® - aye + + - 2 
pe ay 
But the integral 0 has no closed form 


solution. Therefore we have to evaluate it numerically. 
Using the same principle as that in [Ref. 17] the average availability over 


a cycle of length t, + R, and hence in the long run is given by 


] 
ty + R, 





AG) = Г > AD de (4.24) 


By taking the derivative of equation (4.19), we obtain 


dar) ty * R 4) - ft A(t) dt 
dt, . с, + R) 





(4.25) 


Setting 
dA (t,) 
— — — - 0 , fort, » 0 then t,' is a candidate for 
d (t,) a time between the end of one preventive 
replacement and the beginning of the next 


and t, is the value that satisfies the following equation 





1 tp" 
Ай“) = ? A(t) dt 4.2 
(% те, (0 (4.26) 


To solve equation (4.24) the IMSL subroutine DQDAG was used to 
compute the integral in equation (4.23) for various values of t. So we obtained 


A(t) for various values of t and the graph of this function is shown in Figure 
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4.3. Then we used the IMSL subroutine DCSINT [Ref. 13] to compute the cubic 
spline interpolant to the set of data points obtained earlier (values of t and the 
corresponding values of A(t) ). Finally we used the IMSL subroutine DCSITG 
[Ref. 13] to evaluate the integral of the cubic spline for various values of t,. A 
brief description of the IMSL subroutines DCSINT and DCSITG are in 
Appendix D and E respectively. The values obtained by using this ’exact’ 
method was compared with those values obtained using the approximation as 


in equation (4.12). These values are tabulated below. It is observed that the 
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Figure 4.3 Plot of the Availability Function for Values of a= 3.0, 
p = 1390, R, = 8 hours, R, = 8 hours. 





both the approximate and exact solutions are the same to four decimal places 
and the approximation gives an accurate solution. This is true for small values 
of R; and as Rẹ is increased from 1.0 hour to 8 hours, both the approximation 
and the exact solutions are the same to two decimal places, as shown in Table 
4.4. However the optimal replacement age are the same for both the cases. As 
R, increases the accuracy of the approximation diminishes. 

Table 4.3 shows the values of availability for a = 3.0, n = 1390, R, = 17 
hours, R, = 8 hours. The optimal replacement age for the two methods and the 


maximum availability are: 


* approximate solution 2203 hours 0.994581 


* exact solution 2208 hours 0.994589 


It is observed that the approximation gives very accurate results. 
Table 4.4 shows the values of availability г а = 3.0, р = 1390, В. = 8 
hours, R, = 8 hours. The optimal replacement age for the two methods and the 


maximum availability are: 


* approximate solution 1099 hours 0.989201 


* exact solution 1108 hours 0.989301 


It is again observed that the approximation gives very accurate results. 
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Table 4.5 shows the values of availability for a = 2.5, p = 1890, R, = 1 
hour, R, = 8 hours. The optimal replacement age for the two methods and the 


maximum availability are: 


* approximate solution 2711 hours 0.995104 


e exact solution 2717 hours 0.995112 


Once again it is observed that the approximation gives very accurate 
results. If the actual shape parameter was 2.5 and if we have estimated it to 
be 3.0, then we would have lost an availability of 0.000523 which can be 
considered negligible. 

Table 4.6 shows the values of availability for a = 3.5, p = 1390, R, = 1 
hour, R, = 8 hours. The optimal replacement age for the two methods and the 


maximum availability are: 


* approximate solution 1935 hours 0.994243 


* exact solution 1940 hours 0.994253 


Once again it is observed that the approximation gives very accurate 
results. If the actual shape parameter was 2.5 and if we have estimated it to 
be 3.5, then we would have lost an availability of 0.000859 which again can be 
considered negligible. This shows that the shape parameter does not effect the 


availability drastically and a close estimate 1s sufficient. 
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Table 4.7 shows the values of availability for a = 3.0, p = 1350, R, = 1 
hour, R, = 8 hours. The optimal replacement age for the two methods and the 


maximum availability are: 


* approximate solution 2140 hours 0.994421 


* exact solution 2145 hours 0.994430 


It is observed that the approximation gives very accurate results. If the 
actual scale parameter was 1390 and if we have estimated it to be 1350, then 
we would have lost an availability of 0.000159 which again can be considered 
negligible. 

Table 4.8 shows the values of availability for a. = 3.0, p = 1450, R; = 1 
hour, R, = 8 hours. The optimal replacement age for the two methods and the 


maximum availability are: 


* approximate solution 2298 hours 0.994805 


e exact solution 2302 hours 0.994812 


It is observed that the approximation gives very accurate results. If the 
actual scale parameter was 1450 and if we have estimated it to be 1350, then 
we would have lost an availability of 0.000382 which again can be considered 
negligible. This shows that the scale parameter does not effect the availability 


drastically and a close estimate is sufficient. 
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TABLE 4.3 COMPARISON OF AVAILABILITY FOR VALUES OF 
a = 3.0, u = 1390, R; = 1 HOUR, Rp = 8 HOURS. 









Replacement Approximate Exact 
Age (hours) Availability Availability 


БЕСІНШІ ЕНБЕГЕН 


| 2600.0 0.994423 0.994437 | 
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TABLE 4.4 COMPARISON OF AVAILABILITY FOR VALUES OF 
o 2 3.0, u = 1390, R, = 8 HOURS, R, = 8 HOURS. 


Replacement Approximate Exact 
Age (hours) Availability Availability 


T e 1 
1700.0 0.989300 
1800.0 0.986252 
2000.0 0.984530 
0.982519 
2400.0 0.980244 


2600.0 0.976857 | 0.977724 
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TABLE 4.5 COMPARISON OF AVAILABILITY FOR VALUES OF 
a = 2.5, p = 1390, R, = 1 HOUR, R, = 8 HOURS. 


Replacement Approximate Exact 
Age (hours) Availability Availability 


[ wm | sem | sem _ 
o | NN NN 
Сао је | еве 
[ms | vem | omen — 


3400.0 0.994907 0.994921 
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TABLE 4.6 COMPARISON OF AVAILABILITY FOR VALUES OF 
a = 3.5, р = 1390, R, = 1 HOUR, R, = 8 HOURS. 















Replacement Approximate Exact 
Age (hours) Availability Availability 
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TABLE 4.7 COMPARISON OF AVAILABILITY FOR VALUES OF 
a = 3.0, np = 1350, R; = 1 HOUR, R, = 8 HOURS. 















Replacement Approximate Exact 
Age (hours) Availability Availability 


Ошо ЕСТЕН ЕСТЕН 
жо раме ЕСЕН 
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TABLE 4.8 COMPARISON OF AVAILABILITY FOR VALUES OF 
a = 3.0, p = 1450, R, = 1 HOUR, R, = 8 HOURS. 


Replacement Approximate Exact 
Age (hours) Availability Availability 


шш је _ 
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V. POLICY III (MINIMAL REPAIR / FAILURE REPLACEMENT / 


PREVENTIVE REPLACEMENT) 


POLICY II assumes that each time a failure occurs it can be repaired and 
the system is restored to an operational state without changing the failure rate 
of the system as a whole. Now we shall be more realistic and consider the 
possibility of a major failure before the planned age replacement time e. and 
that failure is rectified by a replacement. (Here again a replacement may be 
an overhaul which is assumed to return the system to a state as good as new). 
So when the system is running, two types of failures are possible: 

* Typel failure, denoted by 1C. This failure is corrected by minimal repair; 
if the Type I failure occurs at age t, the failure rate just after correction 
is h(t). This type of failure is usually repaired at sea by the ship’s staff 
using the support elements onboard. 

• Туре П failure denoted by 2C. This failure is remedied by effective system 
replacement; the Type II failure is followed by overhaul, after which the 
failure rate is h(0). This type of failure is beyond the ship’s staff 
capability either due to lack of expertise or unavailability of the required 
support elements, and the ship has to return to base to effect repair 
either by major overhaul or replacement. 

The failures are detected immediately. 


It is assumed that the costs and downtime for simple preventive 


maintenance actions are negligible. 
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First let us assume the time to some type of failure, T, has an IFR 
distribution F(t) with a failure rate denoted by h(t); 0 < t is the system 
operating time since last overhaul. Let Y be the time of a type 2C failure 
having the distribution function G(t). 

Let p, be the probability of Type 1C failure 

and p, be the probability of Type 2C failure 

(p, and p, could dependent on time) ; p, + p, = 1.0 

Now let N(t) denote the number of Type 1C failures that occurs in time 
t, where t is measured from a moment of replacement or overhaul. By the 
assumption, the expected number of Type 1C failures in time t, m(t) 2 E[N(t)], 
satisfies a simple differential equation obtained as follows: 
m(t+dt) = m(t){1 - pıh(t)dt} + {m(t)+1}p,h(t)dt + o(dt) (5.1) 
where o(dt) represents higher order terms which are = 0 
rearranging the terms we get 


dm (t) 





= p,h(t) 
dt (5.2) 


and therefore 


м0 = Py |, Ко 4 (5.3) 
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provided minimal repair correction times are assumed to be negligible enough 
so that the expected downtime from that source is close to R,.-m(t), where R,c 
is the mean time of a minimal repair. 

Let T be the time to a failure (either Type 1C or 2C) with the failure rate 
h(t) 

Then the failure rate of Type I failure = p,h(t) at age t 


and the failure rate of Type II failure = p,h(t) at age t. 


Since 
п-в] = е 9 (5.4) 
and 
1 -од) -е е“ 
ü (е A(t) ty 
1 - GO ] [1 - F( P (5.5) 


From equation (5.5), it can be seen that if p, = 1, then the distribution of 
any type of failure is the same as the distribution of a Type 2C failure and this 
implies that each failure is a major failure and it requires a replacement / 
overhaul action and this is the case as in Policy I (Age Replacement) as 
mentioned in Chapter III. 

For p, = 0, this implies that the probability of a major failure, i.e Type 2C 


is 0 and each failure is a minimal failure and the failure is removed by 
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minimal repair and this is the case as in Policy II (Minimal Repair and 
Preventive Replacement) as mentioned in Chapter IV. 


For the case when 0 < p, < 1, then a failure could be either 


° Type ІС 


° Type 2C 


This is Policy III, which is the general case and will be discussed in this 
chapter. 

This idea of two types of failure has been studied by Beichelt and Fisher 
[Ref. 18]. They derived the Reliability functions for calculating the expected 
long run cost rate for a generalised age-replacement policy. They assumed that 
maintenance actions take only negligible times which in reality is not true, 
especially when a major failure occurs at sea. Besides that, in the military 
environment we are often interested in the availability of the system as our 
measure of effectiveness so in this policy we will expand the cost rate model, 
but also, and more importantly, formulate a model to maximize the availability 


of the system subjected to two types of failures. 


A. COST MODEL 
The age t, at which the operating system is replaced or overhauled 


depends on the following factors: 
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* Failure distribution (IFR). 
* Costs of minimal repair, preventive replacement and failure replacement. 
* Downtime of preventve replacement and failure replacement. 
* Measure of effectiveness: 

- Minimize costs. 

- Maximize availability. 

It is assumed that the downtime of a minimal repair is negligible. 

Let C,. be the cost of a minimal repair. 

С» be the cost of a failure replacement. 

ор be the cost of a preventive replacement. 


c be the mean time of a minimal repair (Type 1C failure). 


O 


t 


c be the mean time of a failure replacement (Type 2C failure). 


Ge 
"Ü 


be the mean time of a preventive replacement. t,. be the planned 
replacement age. 
N(t,) be the number of minimal repairs in interval (0,t,). 
y be the observed time from system replacement until the next 
Type 2C failure. 
N(y) be the number of minimal repairs in the period (0,y), where y is 
the time until Type 2C failure. 


Ric > Rip > Rip and Co > Cop > Ci 


A replacement takes place either at time t, or when there is a Type 2C 


failure. 
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1. Case l 


Y b 
ZR 
[— a a x tH —— 
0 t, 
x----minimal repair 
2. Case 2 
y < t 


P 


Е---К = 
P ВЕН че | — . 
0 y 


x----minimal repair 


A cycle is completed each time a replacement takes place and the 
costs C incurred in a cycle is given by the total costs of minimal repairs and 
the cost of either a preventive replacement or a failure replacement. This is 
given by the following: 

4 | Се [ N(t,) | Y >t, ] + Co w.p. 1 - G(y) 
ОС [Г ну) | AE + С. w.p. G(y) (5.6) 

In order to find the total costs of minimal repairs, we ought to know 
the expected number of minimal repairs and this is given below for both the 


two cases: 
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For case 1, Y » t, 


E[ NG.) | Y » t) - f* p, ht) dr (5.7) 


For case 2, Y < t, 
Let N(y) =n, then 


EI NO) | Y<] = Y, n PI[N(y)7n | Y<t,] 


п=0 (5.8) 
m: P[N(y)=n , Y st,) 
» > Р[Ү = 1] 
thus 
EL NO) | Y «1,1 - —L— f* tf p, &(9 ds 1 g() dr (5.9) 
P G(t,) o Jo ^! 
Then 


FIC] = K(C,c E(NG) | Y * £D * C, II1-G()] * 
(С. EINO) | Y s 5) + Се IGE) 


The length of the cycle, L, depends on the time of a Type 2C failure 


(5.10) 


and is given by: 


( Е. + t, Y ЭА, 
L ={ (5.11) 
{ RetY Ү 


IA 
ct 


Then the expected length of the cycle, E[L], is 


EL] - Rl -G(t)] * RAGE) + MCE) (5.12) 


where M(t,) is the mean age of the system at replacement / overhaul and is 


given by 
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Mt) - f*ts dt + «I1 - GG) ) (5.13) 


which can be simplified to 





ма) - I [1 - СО ] dt (5.14) 
From equation (3.1), 
E[C] 
Cie) a= 
E[L] 


Then the optimal value of Е is the value that minimizes the cost rate 
function C(t,) and this can be found by graphical or numerical analysis. [Ref. 
16] highlights that an finite optimal interval t, exists if h(t) is monotone 
increasing function and C > 0 where C is given by 


Cic 
DN c 
(Cac - С, - Cic) 


3. Weibull Example 
For a time to failure, T, following a Weibull distribution with shape 
parameter a and scale parameter A the expected total cost in a cycle, E[C], is 


given as follows: 
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E[C] - e? (Сер, Ск)“ ч Ср - Се) + Сс 


ляк (5.15) 
Ссрра 224 J: 12% - іе РМ) & 


t ДЕ а 
[г-к 


The integral has no closed form 


solution and the IMSL subroutine DQDAG was used to solve it numerically. 


The length of the cycle is given by 
f РАМ, _ е РЕ ip g PAAD" (5.16) 
ЕЁ = Юре EP) + Ral e] e foe dt 
4. Numerical Illustration 


Now let us consider an example with these data: assuming that the 


downtime for a minimal repair is negligible: 


Cop = $ 25,000 
Сс = $ 37,500 
Cic =$ 1,000 
ЊЕ = 8 hours 


К.с = 16 hours 
Ric = 1 hour 

a = 3.0 

А = 1/1390 hours 
р, = 0.6 


р, = 0.4 
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Figure 5.1 shows the plot of the expected cost rate function. It is 
observed that the curve is fairly flat near the optimal point thus giving 
flexibility to the decision maker. The optimal replacement age is t, = 1888.64 
hours at a cost of $22.03 per hour. 

In a real situation p, is a function of age and it usually increases 
with age i.e. as the system ages the probability of a major failure approaches 
1.0. For our case we have taken the probability to be constant and this again 


is to simplify our computation. 


ГА 
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Figure 5.1 Plot of the Cost Rate Function. 
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TABLE 5.1 VALUES OF COST/HOUR FOR VARIOUS 
REPLACEMENT AGE. 


Replacement age Cost 
t, (hours) ($ / hour) 


тшш men | 
1700.0 22.15 
(mm ae 
= 
аа јав 
a — 
me јав 


2400.0 22.42 


2500.0 
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B. AVAILABILITY MODEL 

Now let us take the (Point) Availability as our measure of effectiveness. 
When the ship is at sea there are usually two types of failures; failures which 
can be rectified onboard, each of which individually does not effect the 
availability of the system drastically, and other failures which are not 
repairable onboard so the ship has to return to port to effect repairs; this 
causes the availability to be more severely degraded. First we shall assume 
thatthe downtime for a minimal repair is negligible. Then using the equation 


(3.8), the availability function is given by 


ма) 
AE RCA (5.17) 
Rp O - G(t)] + Re [GG] + MG) 
where M(t,) is the mean operational time in a cycle and it is given by 
M(t) = ү ? [1 - G()) ] dt (5.18) 


and this is similar to the availability function in Policy I in Chapter III. 
Again the availability function in equation (5.17) is an approximation 

because we do not take into account the downtime during a minimal repair, 

which in reality is not faithful to reality, but should be reasonably accurate. 
If we do take the downtime of the minimal repair into account we can 


obtain an exact solution as follows: 
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Using the same principle as in equation (4.15), 
Let 
Y = Time of a Major Failure 
{1 if the system is up at time t 
X(t)= | 
(0 if the system is down at time t 
B = P{Failure is Type 1С) 
а (б9р ео 
a(t)= P{Y >t, X(t) = 0} 
a,(t) + a(t) = P{Y > t} = 1 - Git) 
{ X(t) if Y>t 
Z(t)= | 
(0 ifYst 
h(t) = the failure rate of Type 1C failure, which is 1 ,(t) 
p(t) = the repair rate of Type 1C failure, and for the Exponential 
distribution, u(t) = Ric? 
Now, 


a, (t+dt) а, (Е) {1 - h(t)dt) * a,(t)u(t)dt + o(dt) (5.19) 


ag(t*dt) - a,(t)(1 - n(t)dt) * a,(t)h(t)dt + o(dt) (5.20) 
Initial conditions: a,(0) 21; a,(0 = 0 
Equation (5.19) becomes 


d (a) 





= - h(t)a,(t) + a,(t)u(t) (5.21) 
dt 
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d (aj) 





= — и(Е)а (Е) + а, (Е)Ъ (+) В (5.22) 
dt 


Since a,(t) = {1 - G(t)} - a,(t), Equation (5.21) becomes 


d (ai) 
“ - h(t)a,(t) + pn(t) ([3 - G(t)] - a,(t)) (5.23) 





dt 
which can be simplified to 


d (a) 





= — (h(t) * w(t) }a,(t) + p(t) {1 - G(t)} (5.24) 
dt 


The equation, 


d(a,) 





+ (h(t) + u(t)} а,(Е) =0 ; with a,(0) = 1 


has the solution 
a,(t) = EXP{- E(t)} 


where 
Е@ = | th) + 165) ds 
Therefore the equation (5.24) has the solution 
a) =e" [1 + ftu) (1.7. GG) e79 ds] (5.25) 


and equation (5.25) can now be written as 
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-J the) +u), -f ht) nias, | 
Р(9 -е — *['ueti-Ge)le '* ^. dt, (5.26) 


and for the Weibull distribution with shape parameter a and scale parameter 
A, and substituting the values for h(t) and p(t) in equation (5.26), the 


probability that the system is up at time t is given by 


2 
Рен -p,G.0* - p,O49* p,049* --- 1 
po = e P RA, L pe mem one He ROO, (Бат 
1 R, Jo 2 
Ip AD" - par) da 
; PO A + PAD t OTR 
f e с Яс j, 
But the integral “0 has no closed-form 


expression in terms of elementary tabulated functions. Therefore we have to 
evaluate it numerically. 


Since the length of the cycle, L is 


{ Кор + t, Y > to 
L ={ 
{ Re+ Y Y<t, 
then using the same principle as that used in [Ref. 16] the average availability 


over a cycle of length L and hence in the long run is given by 


t, 
ад (5.28) 
Expected length of a Cycle 


But the Expected length of a cycle is given by equation (5.12). Therefore 


A(t) - 


the Availability, A(t,), is given by 
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i p a,(t) dt 
O RR (5.29) 
Rl - G(t)) * R,L(G(t)) + i^ 1 - G(t) dt 

To study the above expression the IMSL subroutine DQDAG was used to 

compute the Availability function a,(t) for various values of t. Then we used 

the IMSL subroutine DCSINT [Ref. 13] to compute cubic spline interpolant to 

the set of data points obtained earlier (values of t and the corresponding values 

of a,(t) ). Finally we used the IMSL subroutine DCSITG [Ref. 13] to evaluate 

the integral of the cubic spline for various values of the replacement age t.. 

From this we can compute the optimal replacement age t, and the maximum 

availability. The values obtained by using this ’exact’ solution was compared 

with the approximation as in equation (5.17). These values are tabulated 

below. It is observed that both the approximate and exact solutions are the 

same to three decimal places, and the approximation gives an accurate 

solution. However this is true for small values of R,. and R,.. As the values of 

Rc and Rc are increased the accuracy of the present approximation method 
degrades. 

Table 5.2 shows the values of availability for & 2 3.0, n 2 1390, Ri, 2 1 

hr, R,. = 16 hrs, R;, = 8 hrs. The optimal replacement age for the two methods 


and the maximum availability are: 
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* approximate solution 1528 hours 0.991715 


* exact solution 1496 hours 0.991567 


If we use the approximate solution as our reference for the replacement action, 
the exact availability is 0.991565, which is extremely close to the exact 
solution, and much more easily obtained. 

Table 5.3 shows the values of availability for a = 3.0, p = 1390, R,, = 8 
hr, R,. = 24 hrs, R,p = 8 hrs. The optimal replacement age for the two methods 


and the maximum availability are: 


* approximate solution 1201 hours 0.989796 


* exact solution 1113 hours 0.989070 


If we use the approximate solution as our reference for the replacement action, 
the exact availability is 0.989015 which is extremely close to the exact solution. 

Table 5.4 shows the values of availability for a = 3.0, u = 13890, R,, = 8 
hr, R,.. = 72 hrs, R,p = 8 hrs. The optimal replacement age for the two methods 


and the maximum availability are: 


* approximate solution 752 hours 0.984143 


e exact solution 736 hours 0.983835 


If we use the approximate solution as our reference for the replacement action, 


the exact availability is 0.983827 which is very close to the exact solution. 
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Table 5.5 shows the values of availability for a = 3.5, p = 1390, Вис = 8 
hr, R,. = 24 hrs, R.p = 8 hrs. The optimal replacement age for the two methods 


and the maximum availability are: 


* approximate solution 1147 hours 0.990119 


* exact solution 1076 hours 0.989519 


If we use the approximate solution as our reference for the replacement action, 
the exact availability is 0.989471 which is again very close to the exact 
solution. 

Table 5.6 shows the values of availability for a = 2.5, p = 1390, R,, = 8 
hr, R,c = 24 hrs, R,p = 8 hrs. The optimal replacement age for the two methods 


and the maximum availability are: 


* approximate solution 1317 hours 0.989493 


* exact solution 1200 hours 0.988594 


If we use the approximate solution as our reference for the replacement action, 
the exact availability is 0.988535 which is once again very close to the exact 
solution. 

From the above results we can conclude that the approximate solution is 
a good approximation for planning the replacement or overhaul actions of a 


system in order to maximize the availability of the system. From our earlier 
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discussion we found that the availability is sensitive to the shape parameter 
a; Figure 5.2 shows the availability (exact solution) for various values of alpha. 

If the shape parameter a. was in fact 3.5 and in our estimation we used 
a = 2.5, then we would have lost an availability of 0.000626, which is 
extremely small. This shows that the parameters need not be estimated very 
accurately to achieve good results. 

Figure 5.2 shows that for low values of alpha such as 2.5, the availability 
function is rather flat and there is more flexibility in determining the 
replacement age that is, the replacement interval at 1200 hours or 1500 hours 
gives about the same availability on the long run. However this is not true for 
higher values of alpha such as 3.5 where the availability function falls quite 
rapidly on both sides of the optimal replacement age t,'. At the availability of 
about 0.9883, the replacement age is 1475 hours and it is insensitive to the 
value of alpha. 

Now we shall also look at the effects of the scale parameter A on the 
availability. Table 5.7 shows the values of availability for a = 3.0, p = 1350, R,, 
= 8 hr, Кос = 24 Бтв, Кр - 8 hrs. The optimal replacement age for the two 


methods and the maximum availability are: 


* approximate solution 1167 hours 0.989496 


* exact solution 1081 hours 0.988750 
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Figure 5.2 The Plot of ‘Exact’ Availability for Various Values of the 
Shape Parameter Alpha. 





If we use the approximate solution as our reference for the replacement 
action, the exact availability is 0.988693 which is very close to the exact 
solution. 

Now if the actual scale parameter was 1390 hours and it was estimated 
to be 1350 hours, then we would have lost an availability of 0.00032 which can 
be considered very small. 

Table 5.8 shows the values of availability for a = 3.0, p = 1450, R,, = 8 
hr, R,. = 24 hrs, R,p = 8 hrs. The optimal replacement age for the two methods 


and the maximum availability are: 
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TABLE 5.2 COMPARISON OF AVAILABILITY FOR VALUES OF 
a = 3.0, a = 1390, R,. = 1 HR, R,. = 16 HRS, Rp = 8 HRS 


Replacement Approximate Exact 
Age (hours) Availability Availability 


| 13200 СЕС 0991592 [o — =“ 
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TABLE 5.3 COMPARISON OF AVAILABILITY FOR VALUES OF 


Replacement Approximate Exact 
Age (hours) Availability Availability 


me |мен | on 


| 1354.0 0.989677 0.988710 
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TABLE 5.4 COMPARISON OF AVAILABILITY FOR VALUES OF 
a = 3.0, u = 1390, R,. = 8 HRS, E.... = 72 HRS, Rp = 8 HRS 


Replacement Approximate Exact 
Age (hours) Availability Availability 


жо Ге онно 
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TABLE 5.5 COMPARISON OF AVAILABILITY FOR VALUES OF 
а = 3.5, u = 1390, R,. = 8 HRS, R,. - 24 HRS, R,, - 8 HRS 

















Approximate Exact 
Age (hours) Availability Availability 


[x 

Cmm E ЕН ЕСЕН 
1208.0 0.989361 

ЕСІЕЕН ЕСЕН ЕСЕН 


| 1280.0 0.989994 0.989160 | 


Replacement 
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TABLE 5.6 COMPARISON OF AVAILABILITY FOR VALUES OF 
a = 2.5, u = 1390, R.. =8 HRS, R,. = 24 HRS, R,, = 8 HRS 


Replacement Approximate Exact 
Age (hours) Availability Availability 


1456.0 0.989433 0.988342 
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TABLE 5.7 COMPARISON OF AVAILABILITY FOR VALUES OF 
a= 3.0, p= 1350, R,. = 8 HRS, R,. = 24 HRS, R,r = 8 HRS 


Replacement Approximate Exact 
Age (hours) Availability Availability 


m [s 


0.989446 0.988523 
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TABLE 5.8 COMPARISON OF AVAILABILITY FOR VALUES OF 
a = 3.0, p = 1450, R,, = 8 HRS, R;, - 24 HRS, R,, - 8 HRS. 













Replacement Approximate Exact 
Age (hours) Availability Availability 


[ww [sem Tran _ 
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* approximate solution 1253 hours 0.990214 


* exact solution 1161 hours 0.989516 


If we use the approximate solution as our reference for the replacement 
action, the exact availability is 0.989464 which is again very close to the exact 
solution. 

Now if the actual scale parameter was 1450 hours and it was estimated 
to be 1350 hours, then we would have lost an availability of 0.000446 which 
once again can be considered very small. This shows that the scale parameter 


does not effect the availability drastically and a close estimate is sufficient. 


C. MAXIMIZE AVAILABILITY SUBJECT TO BUDGET CONSTRAINT 

The availability and cost are very important measures of effectiveness. 
We would like to have as many resources as possible to maximize the 
availability of a system, however in reality we are often limited by budget 
constraints. So we would like to achieve cost effectiveness, that is we would 
like to 

maximize Effectiveness Level 

subject to Budget < B 

The effectiveness level which is commonly used is the availability. So in 
our case we would like to 

Maximize A(t,) 


Subject to СС) < B. 
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We have solved this problem graphically. Figure 5.3 shows the plots of the 
availability and the cost rate functions as formulated in Section B. It 1s 
observed that the availability is maximum at 0.991567 at an optimal 
replacement interval ES of 1496 hours while the cost is minimum at 
$22.03/hour at an optimal replacement interval t, of 1889 hours. 

Now if we have a budget of not exceeding $22.25/hour, then the maximum 
availability that can be obtained is 0.9915 with an optimal replacement 


interval t, of 1650 hours. 


— Cost /} hour 


= Avai lability 


- — red! 


1200 1500 1800 2100 2400 2700 3000 


Time Chours) 


Figure 5.3 Plot Showing the Availability and the Cost Rate 
Function. 





The other measure of effectiveness that is also used is the mission 


reliability and this again will give another optimal replacement age, so now we 
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have a multiple independent conflicting criterion and it is up to the Decision 
Maker to decide which of the measure of effectiveness is vital and fits the 


scenario very well. 
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VI. CONCLUSION AND RECOMMENDATIONS 


The fighting effectiveness and operational readiness of a ship depends 
largely on the operational availability of her equipment and systems. If we are 
not constrained by the budget then we can expend the resources necessary to 
achieve the desired availability, such as incorporating more redundancies or 
carrying out "premature" replacements and overhauls. However in reality this 
is not the case, and we are always limited by the available budget, and we 
would like to maximize system availability subject to budget constraints. 

At sea, equipment and systems are exposed to various environment and 
unfavourable conditions. As such, they are subject to stochastic failure and 
deterioration. However, with timely maintenance actions as discussed in the 
various policies in the thesis, we can minimize catastrophic and unexpected 
failures, enabling us to achieve the desired measures of effectiveness. 
Therefore, based on the optimal maintenance policies, we can carry out 
replacement actions or complete overhaul of the equipments and systems at 
the base during the stand-off periods so that when the ship is at sea we can 
minimize loss of availability due to failures and maximize our successes in the 
operational missions. 

Based on the policies, we can also carry adequate spares. This is 


particularly important for long missions or when the need to be prepositioned 
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in forward operating areas. On the other hand, although we can carry spares 
onboard, sometimes other support elements are not available to rectify the 
defects. 

In our analysis, we have assumed that we know the failure distribution. 
For our case we concentrated on the Weibull distribution. The two important 
parameters which are usually estimated are the shape parameter, alpha, and 
the scale parameter, lambda. From the analysis it is observed that these 
parameters need not be estimated very accurately. A slight variation in the 
values of the parameters do not affect the long run availability and cost rate 
functions drastically. Of the two, the shape parameter needs to be estimated 
more accurately. 

In the formulation of the availability functions, we simplified the 
computation by assuming the downtime for a minimal repair to be negligible. 
When compared with the exact solution, taking into account the downtime for 
a minimal repair, the results obtained by the approximation method gave 
extremely accurate results. Many of the functions and integrals that were 
formulated did not have closed form expressions in terms of elementary 
tabulated functions. However IMSL subroutines were available in the Math 
library at the main frame at the Naval Postgraduate School and these 
subroutines expeditiously computed the integrals very accurately. 

The expected downtimes, as taken in our analysis, are practical figures 


assuming all the support elements are readily avalilable when required. That 
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is the reason for the high values of the availability. In reality we often have to 
wait for spare par s and support elements, which sometimes have long lead 
times. 

In the thesis we also studied the effects of simple preventive maintenance 
actions on the measures of effectiveness. It is observed that simple preventive 
maintenance actions do not restore the system to a condition "as good as new" 
but the maintenance actions can enhance or improve the reliability of the 
system by a certain factor which decreases as the system ages. However, a 
methodology for chacterizing the effective age reduction remains to be 
developed. In our analysis we have assumed that the system improves by a 
certain factor on completion of each preventive maintenance action. 

Most planned maintenance systems usually adopt maintenance efforts 
based on calendar time (weeks, months) or running hours of systems or 
equipment, but from the analysis it is observed that for systems that have a 
"wear-out' life distribution we shall have to successively resort to decreasing 
maintenance intervals if we are going to maintain the systems above some 
minimum reliability level. 

This subject can be expanded further by future research. The following 
areas are recommended: 

* Carry out similar analysis, especially for the availability function for the 


three policies when the underlying life distribution F comes from a family 
of Gamma Distribution. 
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е Incorporate imperfect repair into the model. 


* À methodology for chacterizing the improvement of the reliability of the 
system on completion of preventive maintenance needs to be investigated. 


* À system may consists of various sub-systems. The maintenance actions 
for the sub-systems has to be coordinated so that instead of just taking 
into account the availability of the sub-systems individually, the whole 
system has to be considered bearing in mind of some dependence on the 
supporting elements associated with the sub-systems. 
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APPENDIX A. COMPUTER PROGRAM AND OUTPUT 


PROGRAM PM 


THIS PROGRAM IS FORMULATED TO COMPUTE THE AVAILABILITY, 
RELIABILITY AND COST FUNCTIONS FOR VARIOUS POLICIES OF 
INTEREST. THE PROGRAM DETERMINES THE OPTIMAL REPLACEMENT 
OR OVERHAUL INTERVALS FOR THE SYSTEM. THE OPTIMAL 
REPLACEMENT INTERVAL DEPENDS ON THE MEASURES OF 
EFFECTIVENESS DESIRED AND THE RELEVANT TYPE OF POLICY 
APPLICABLE TO THE SYSTEM UNDER STUDY. THE PROGRAM 
UTILIZES THE IMSL SUBROUTINES FROM THE MATH LIBRARY 
AVAILABLE AT THE MAIN FRAME AT THE NAVAL POST GRADUATE 
SCHOOL. IT IS ASSUMED THAT THE TIME TO FAILURE FOLLOWS 
THE WEIBULL DISTRIBUTION WITH SHAPE PARAMETER ALPHA AND 
SCALE PARAMETER LAMBDA. 


Фе Фо е оо Фо ООС 


PARAMETER STATEMENT 
INTEGER NDATA,NTINT 


PARAMETER (NDATA = 81, NTINT = 20) 


C VARIABLES DECLARATION 


REAL*8 A, B, F, RESULT, ERRABS, ERRREL, ERREST, ALPHA, 
LAMBDA, MU, C2C, C1C, C2P, R2P, R2C, R1P, R1C, 

Al, AA, XDATA(NDATA), YDATA(NDATA), 
BREAK(NDATA), CSCOEF(4,NDATA), AVAIL, AVAIL1, 

P1, P2, SOLN, SF, CDF, LENGTH, G, H, ANSWER, 
DCSITG, COST, COSMIN, ANS, D, REL, E, EE, FF, 

TD2P, TDT1, TDT2, SUMD1P, SUMD2C, STIME, TL, 
CDFTL, T(NTINT), REL, R, TDT, UPTIME, TDIP, 

SUM, SUMI1, X, XX, LL, UL, STEP 


Со Ro Ro Ro Ro Ro Ro Ro 


INTEGER IRULE,POL,P,SPMREQ,NN 


C VARIABLES DEFINITION 
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(C2 QC) 021022 COSE COL CO CO" C2 CO" COT ОЕЕО оо Ое ооо ообо 


Arm... LOWER LIMIT OF INTEGRATION USED IN 
ARGUMENT OF THE IMSL SUBROUTINE DQDAG AND 
IMSL FUNCTION DCSITG 

pos UPPER LIMIT OF INTEGRATION USED IN 
ARGUMENT OF THE IMSL SUBROUTINE DQDAG AND 
IMSL FUNCTION DCSITG 

E,EE,F,FF,G,H..FUNCTIONS TO BE INTEGRATED 

RESULT......... ESTIMATE OF THE INTEGRAL FROM A TO B OF 
THE FUNCTIONS E,EE,F,FF 

ERRABS......... ABSOLUTE ACCURACY DESIRED AS THE INPUT 
ARGUMENT OF THE IMSL SUBROUTINES DQDAG 

ERRREL......... RELATIVE ACCURACY DESIRED AS THE INPUT 
ARGUMENT OF THE IMSL SUBROUTINES DQDAG 

ЕО... CHOICE OF QUADRATURE RULE (SEE APPENDIX 


ALPHA.......... SHAPE PARAMETER OF WEIBULL DISTRIBUTION 
LAMBDA......... SCALE PARAMETER OF WEIBULL DISTRIBUTION 


NI ........... RECIPROCAL OF LAMBDA 

ЕЛЕНЕ... EXPECTED COST OF A FAILURE REPLACEMENT 

ВАС... EXPECTED COST OF A MINIMAL REPAIR 

ӘДЕТ... EXPECTED COST OF A PREVENTIVE REPLACEMENT 

ECC, eoo EXPECTED DOWNTIME OF A FAILURE 
REPLACEMENT 

ВНЕС... сте... EXPECTED DOWNTIME OF A MINIMAL REPAIR 

RAR... EXPECTED DOWNTIME OF A PREVENTIVE 
REPLACEMENT 

ЛЕТТІ... COMPUTE EXPONENTIAL EXPRESSION (POLICY 
ID 

ДАЛАНЫ... AVAILABILITY AT TIME T (POLICY II) 


NDATA.......... NUMBER OF DATA POINTS FOR COMPUTING CUBIC 
SPLINE INTERPOLANT. INPUT ARGUMENT FOR 
IMSL SUBROUTINE DCSINT 

XDATA( )...... ARRAY OF LENGTH NDATA CONTAINING THE DATA 
POINTS ABSCISSAS. INPUT ARGUMENT FOR IMSL 
SUBROUTINE DCSINT 

YDATAC )..... ARRAY OF LENGTH NDATA CONTAINING THE DATA 
POINTS ORDINATES. INPUT ARGUMENT FOR IMSL 
SUBROUTINE DCSINT 

BREAK( )..... ARRAY OF LENGTH NDATA CONTAINING THE 
BREAKPOINTS FOR THE PIECEWISE CUBIC 
REPRESENTATION. OUTPUT ARGUMENT OF THE 
IMSL SUBROUTINE DCSINT 

CSCOEF......... MATRIX OF SIZE 4 BY NDATA CONTAINING THE 
LOCAL COEFFICIENTS OF THE CUBIC PIECES. 
OUTPUT ARGUMENT OF THE IMSL SUBROUTINE 
DCSINT 
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ооо уос о оссо о оо зо еее се 


С) С) 


AVAIL.......... APPROXIMATE AVAILABILITY 
ANSWER......... EXACT AVAILABILITY 


P1..... s PROBABILITY OF TYPE I (MINIMAL) FAILURE 
PZA PROBABILITY OF TYPE II (MAJOR) FAILURE 
Sb... mm SURVIVAL FUNCTION 

СВЕ. иа DISTRIBUTION FUNCTION 


LENGTH......... LENGTH OF A CYCLE 
COST, COSMIN...COST RATE FUNCTION 
¡O em MISSION DURATION 


TD2P,TDT2..... TOTAL DOWNTIME OF PREVENTIVE REPLACEMENT 

TD1P,TDT1..... TOTAL DOWNTIME OF PREVENTIVE MAINTENANCE 

SUMDIP ........ TOTAL DOWNTIME WHEN TIME OF FAILURE IS 
AFTER PREVENTIVE REPLACEMENT 

SUMDZC ........ TOTAL DOWNTIME WHEN TIME OF FAILURE IS 
BEFORE PREVENTIVE REPLACEMENT 

STIME ......... TOTAL DOWNTIME IN A CYCLE 

UPTIME TOTAL UPTIME IN A CYCLE 


DE a es LOWER LIMIT OF REPLACEMENT INTERVAL 
UB. сс: UPPER LIMIT OF REPLACEMENT INTERVAL 
STEP... аи STEP SIZE OR INCREMENT 
РО а POLICY YO BE EVALUATED 


SPMREQ......... SIMPLE PREVENTIVE MAINTENANCE REQUIREMENT 


EXTERNAL SUBROUTINES AND FUNCTIONS 


EXTERNAL E,EE,F,FF,G,H,DQDAG,DCSINT,DCSITG 


COMMON BLOCKS 


COMMON /FT/ R1C,B 

COMMON /GT/ P2 

COMMON /IT/ P1 

COMMON /HT/ LAMBDA,ALPHA 


A = 0.0D0 
ERRABS = 0.0D0 
ERRREL = 0.00001D0 
IRULE = 2 

PREL = 1.000 
SUMD1P = 0.0D0 
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SUMD2C = 0.0D0 
STIME = 0.0D0 
SUM =0.0D0 


Prompt the user for the policy to be evaluated 


С 
C J POLICY I ------- AGE REPLACEMENT 

C POLICY II ------ AGE REPLACEMENT WITH MINIMAL REPAIR 
C POLICY II ----- AGE REPLACEMENT WITH TWO TYPES OF 

C FAILURES. TYPE I MINIMAL FAILURE, TYPE 

C II MAJOR FAILURE 


PRINT*, ENTER THE NUMBER OF THE POLICY TO BE EVALUATED’ 


READ*,POL 
PRINT*,’POLICY = ’,POL 
GOTO (100,200,300),POL 


100 CONTINUE 


PRINT*/DO YOU WANT TO INCLUDE SIMPLE PREVENTIVE 
& MAINTENANCE-ENTER 1 FOR NO AND 2 FOR YES’ 
READ*,SPMREQ 

GOTO (125,150),SPMREQ 


МЕ, С а 
C POLICY I (AGE REPLACEMENT-NO SIMPLE PREVENTIVE 

C MAINTENANCE) 

РЕКЕ ОЛИО D l u. e mn 


125 CONTINUE 
PRINT*,ENTER THE ESTIMATED VALUES FOR THE SHAPE 
& PARAMETER ALPHA AND THE SCALE PARAMETER MU ' 
READ*,ALPHA,MU 
PRINT*,,ALPHA =’,ALPHA,’MU = МО 


PRINT*, ENTER ESTIMATED VALUES FOR THE EXPECTED COST OF 
& PREVENTIVE REPLACEMENT FOLLOWED BY THE EXPECTED 

& COST OF FAILURE REPLACEMENT 

READ*,C2P,C2C 

PRINT* EXPECTED COST OF A PREVENTIVE REPLACEMENT = 

& ',C2P 

& » EXPECTED COST OF A FAILURE REPLACEMENT = 

& „С2С 


PRINT*, ENTER VALUES FOR THE EXPECTED DOWNTIME OF 
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& PREVENTIVE REPLACEMENT FOLLOWED BY EXPECTED 

& DOWNTIME OF A FAILURE REPLACEMENT 

READ*,R2P,R2C 

PRINT*, EXPECTED DOWNTIME OF A PREVENTIVE REPLACEMENT = 
& ’В2Р 

&  , EXPECTED DOWNTIME OF A FAILURE REPLACEMENT = 

& ',R2C 


PRINT*, ENTER THE VALUE OF MISSION DURATION TIME IN 
& HOURS' 

READ*,D 

PRINT*, MISSION DURATION TIME -', D,'HOURS 


PRINT*,ENTER THE VALUE FOR LOWER LIMIT, UPPER LIMIT, 
& STEP SIZE ’ 

READ*,LL,UL,STEP 

PRINT*; LOWER LIMIT = ’, LL, UPPER LIMIT = ’,UL,STEP = 

& "STEP 


WRITE(30,101)ALPHA,MU,C2P,C2C,R2P,R2C,D 
101 FORMAT(15X'POLICY,3X'T,3X (AGE REPLACEMENTY 
&  /,15X,30(-’)/, 
& 10X’ALPHA = ?,ЕЗ.1,17Х, МО = ’F6.1/, 
& 10X,C2P = $,F7.1,1%X, 'C2C = $,F7.1/, 
& 10X,R2P = ’,F4.1’ HOURS’,11X/R2C = ';F41; 
& HOURS 
& ’/,10X, MISSION DURATION = ’F4.1’ HOURS’//, 
& 5X/REPLACEMENT',5X,COST RATE’,5X,’AVAILABILITY 5X, 
& MISSION RELIABILITY /, 
& 5X, AGE (HOURS),5X,($/HOURY,10X,(AY, 17X,(R)/, 
& 5X,110-)|  5X,8(-) 6Х,12(-), 5Х,19(-)) 


LAMBDA = 1.0D0/MU 

DO 11 I= LL,UL,STEP 
B-DBLE(I) 
CALL DQDAG (E,A,B,ERRABS,ERRREL,IRULE,RESULT,ERREST) 
SF=EXP(-((LAMBDA*B)**ALPHA)) 
CDF=1.0D0-SF 
COST=((C2P*SF)+(C2C*CDF))/((R2P*SF)+(R2C*CDF)+RESULT) 
AVAIL=RESULT/((R2P*SF)+(R2C*CDF)+RESULT) 
REL= DEXP(-((LAMBDA*(B+D))**ALPHA)) / 

& DEXP(-((LAMBDA*B)**ALPHA)) 


WRITE(30,102)B,COST,AVAIL,REL 
102 FORMAT(5X,F7.1,9X,F8.3,8X,F8.6,12X,F8.6) 
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11 CONTINUE 


GOTO 99 
IM LLL LL nc ii 
C POLICY I (AGE REPLACEMENT - WITH SIMPLE PREVENTIVE 
C MAINTENANCE) 

В м... aaee 


150 CONTINUE 


PRINT*, ENTER THE ESTIMATED VALUES FOR THE SHAPE 

& PARAMETER ALPHA AND THE SCALE PARAMETER MU’ 
READ*,ALPHA,MU 

PRINT*'ALPHA z', ALPHA,'MU - ', MU 


PRINT*'ENTER VALUES FOR THE EXPECTED DOWNTIME OF 

& PREVENTIVE REPLACEMENT FOLLOWED BY EXPECTED 

& DOWNTIME OF A FAILURE REPLACEMENT 

READ*,R2P,R2C 

PRINT* EXPECTED DOWNTIME OF A PREVENTIVE REPLACEMENT = 
& ',R2P 

&  ,'EXPECTED DOWNTIME OF A FAILURE REPLACEMENT = 

& ',R2C 


PRINT*,ENTER THE EXPECTED DOWNTIME FOR SIMPLE 

& PREVENTIVE MAINTENANCE AND THE IMPROVEMENT 

& FACTOR ON COMPLETION OF THE MAINTENANCE’ 
READ*,R1P,R 

PRINT*, THE DOWNTIME FOR SIMPLE PREVENTIVE MAINTENANCE 
& = RIP 

& , THE IMPROVEMENT FACTOR = R 


PRINT*/ENTER THE VALUE FOR LOWER LIMIT, UPPER LIMIT, 
& STEP SIZE’ 

READ*,LL,UL,STEP 

PRINT*, LOWER LIMIT = ’,LL, UPPER LIMIT =’,UL,STEP = 

& "STEP 


WRITE(30,103)ALPHA,MU,R2P,R2C,R1P,R 
108 FORMAT(SX,POLICY',2X, ''2X, (AGE REPLACEMENT WITH 
& SIMPLE’,1X, PREVENTIVE MAINTENANCE) ,8X, 
« 63 (-), //, 
& 10X,ALPHA = ’,F3.1,17X,MU = ”,F6.1/, 
& 10X,R2P =’,F4.1,” HOURS,11X,R2C = ”,F4.1, 
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& HOURS’,/,10X,R1P =’,F4.1, HOUR’,12X,R = 

& ’,F4.1/, 

& 5X,’NO. OF INTERVAL’, 5X, REPLACEMENT AGE (HOURS),5X, 
& “AVAILABILITY ’,/,5X,15(’-’),5X,23(-”),5X,12(’-’)) 


LAMBDA = 1.0D0/MU 

T(1) = (1.0D0/LAMBDA)*((-(LOG(0.75D0)))** 

& (1.0D0/ALPHA)) 

CDF =1.0D0 - EXP(-((LAMBDA*T(1))**ALPHA)) 
SF = 1.0D0-CDF 

B = T(1) 


DO 10 I = 1,NTINT 
IF (1.GT. 1) THEN 
T(I) =T(I-1)+(((1.0D0-R)**(I-1))*T(1)) 
END IF 
10 CONTINUE 


DO 20 N = LL,UL,STEP 
IF (N .EQ. 1) THEN 
TDT = (R2P*SF)+(R2C*CDF) 
А = 0.000 
CALL DQDAG(FF,A,B,ERRABS,ERRREL,IRULE,RESULT, 
& ERREST) 
UPTIME = (ALPHA*(LAMBDA**ALPHA)*RESULT)+(SF*T(1)) 
ELSE 
DO 30 J = 1,N 
IF (J .EQ. 1) THEN 
TL = 0.0D0 
ELSE 
TL = T(1)-(T(J)-T(J-1)) 
END IF 
CDFTL = 1-EXP(-((LAMBDA*TL)**ALPHA)) 
REL = 1-(CDF-CDFTL) 
PREL = PREL*REL 
30 CONTINUE 


TD1P = (N-1)*R1P*PREL 
TD2P = R2P*PREL 
TDT1 = TD1P+TD2P 
PREL = 1.0D0 


DO 40 K= 2, N 
DO 50 L = 1,К-1 
IF (L .EQ. 1) THEN 
TL = 0.0D0 


114 


ELSE 
TL = T(1)-(T(L)-T(L-1)) 
END IF 
CDFTL = 1-EXP(-((LAMBDA*TL)**ALPHA)) 
REL = 1-(CDF-CDFTL) 
PREL = PREL*REL 
50 CONTINUE 


TL = T(1)-(T(K)-T(K-1)) 
CDFTL = 1-EXP(-((LAMBDA*TL)**ALPHA)) 
FT = CDF-CDFTL 
SUMD1P = SUMD1P+((K-1)*R1P*PREL*FT) 
SUMD2C = SUMD2C+(R2C*PREL*FT) 
PREL = 1.0D0 

40 CONTINUE 


TDT2 = SUMD1P+SUMD2C+(R2C*CDF) 
TDT = TDT1+TDT2 
PREL = 1.0D0 
DO 60 M =2,N 
DO 70 P = 1,M-1 
IF (P .EQ. 1) THEN 
TL = 0.0D0 
ELSE 
TL = T(1)-(T(P)-T(P-1)) 
END IF 
CDFTL = 1-EXP(-((LAMBDA*TL)**ALPHA)) 
REL = 1-(CDF-CDFTL) 
PREL = PREL*REL 
70 CONTINUE 


A = T(1)-(T(M)-T(M-1)) 
CALL DQDAG(FF,A,B,ERRABS,ERRREL,IRULE,RESULT, 
& ERREST) 
ANSWER = ALPHA*(LAMBDA**ALPHA)*RESULT 
TL = T(1)-(T(M)-T(M-1)) 
CDFTL = 1-EXP(-((LAMBDA*TL)**ALPHA)) 
REL = 1-(CDF-CDFTL) 
ANSWER = ANSWER+(REL*(T(M)-T(M-1))) 
STIME = STIME +(PREL*ANSWER) 
PREL = 1.0D0 
60 CONTINUE 


A = 0.0D0 


CALL DQDAG (FF,A,B,ERRABS,ERRREL,IRULE, 
& RESULT,ERREST) 
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ANSWER = (ALPHA*(LAMBDA**ALPHA)*RESULT) + 
$ (SF*T(1)) 

UPTIME = STIME+ ANSWER 
END IF 


AVAIL = UPTIMEAUPTIME+TDT) 
WRITE(80,104)N, TIN), AVAIL 
20 CONTINUE 


104 FORMAT(10X,13,19X,F6.1,17X,F8.6) 
GOTO 99 


C POLICY II (AGE REPLACEMENT - WITH MINIMAL REPAIR) 


200 CONTINUE 
PRINT*,/ENTER THE ESTIMATED VALUES FOR THE SHAPE 
& PARAMETER ALPHA AND THE SCALE PARAMETER MU ' 
READ*,ALPHA,MU 
PRINT*,’ALPHA =’,ALPHA,MU = ’,MU 


PRINT*,ENTER ESTIMATED VALUES FOR THE EXPECTED COST OF 
& PREVENTIVE REPLACEMENT FOLLOWED BY THE EXPECTED 
& COST OF MINIMAL REPAIR’ 

READ*,C2P,C1C 

PRINT*/EXPECTED COST OF A PREVENTIVE REPLACEMENT = 
& *,С2Р 
& ,» EXPECTED COST OF A MINIMAL REPAIR = 
& „СЪС 


PRINT*,,ENTER VALUES FOR THE EXPECTED DOWNTIME OF 

& PREVENTIVE REPLACEMENT FOLLOWED BY EXPECTED 

& DOWNTIME OF A MINIMAL REPAIR” 

READ*,R2P,R1C 

PRINT*, EXPECTED DOWNTIME OF A PREVENTIVE REPLACEMENT = 
& ”R2P 

&  ,'EXPECTED DOWNTIME OF A MINIMAL REPAIR = 

& „ВАС 


PRINT* ENTER THE VALUE FOR LOWER LIMIT, UPPER LIMIT, 
& STEP SIZE’ 

READ*,LL,UL,STEP 

PRINT*, LOWER LIMIT = ’,LL, UPPER LIMIT =’,UL, STEP = 

& STEP 


116 


LAMBDA = 1.0D0/MU 


PRINT*,DO YOU WANT TO INCLUDE SIMPLE PREVENTIVE 
& MAINTENANCE-TYPE 1 FOR NO AND 2 FOR YES' 
READ*,SPMREQ 


GOTO (225,250), SPMREQ 


C POLICY II (AGE REPLACEMENT - WITH MINIMAL REPAIR WITHOUT 
С SIMPLE PREVENTIVE MAINTENANCE) 


225 CONTINUE 


WRITE(30,105)ALPHA,MU,C2P,C1C,R2P,R1C 
105 FORMAT(9X,POLICY',2X,'II,2X, (AGE REPLACEMENT WITH 
& MINIMAL'1X,REPAIR)/,9X,50(-)/, 
& 10X,ALPHA = ’,F3.1,17X,MU = ’,F6.1/, 
& 10X,C2P = $',F7.1,12X,/C1C = $,F7.1/, 
& lOX,R2P =’,F4.1, HOURS’11X,R1C = ’,F4.1, 
& HOURS’,//,5X, REPLACEMENT ’,6X,’COST’,10X, APPROXIMATP’, 
& 8X EXACT'/, 5X,AGE HOURS) 5X,$/HOUR”,9X, ” 
& AVAILABILITY’, 4X, “AVAILABILITY /,5X,11(’-’), 
& 5X,7(-), 8X,12(-),4X,12(-))) 


DO 12 I = 0,4000,50 
B=DBLE(I) 
CALL DQDAG (EE,A,B,ERRABS,ERRREL,IRULE,RESULT,ERREST) 
Al =DEXP(-(((LAMBDA*B)**ALPHA)+(B/R1C))) 
AA = A1+(RESULT*(1.0DO/R1C) ) 
WRITE(13,*)B,AA 
12 CONTINUE 


REWIND (13) 
DO 22 J = 1, NDATA 
READ(13,*)XDATA(J), YDATA(J) 
22 CONTINUE 
CALL DCSINT(NDATA,XDATA, YDATA,BREAK,CSCOEF) 


С 
С Calculate the integral of the spline 
C approximation. 
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DO 32 K = LL,UL,STEP 

B = DBLE(K) 
NINTV = NDATA- 1 
ANSWER = DCSITG(A,B,NINTV,BREAK,CSCOEF) 
ANSWER = ANSWER(B+R2P) 
AVAIL = (B-(R1C*((LAMBDA*B)**ALPHA)))(B+R2P) 
AVAIL1 = BAB+R2P+(R1C*((LAMBDA*B)**ALPHA))) 
COST = ((C1C*((LAMBDA*B)**ALPHA))+C2P)/(B+R2P) 
WRITE (30,106) B,COST,AVAIL,ANSWER 

32 CONTINUE 

106 FORMAT(7X,F6.1,6X,F9.4,10X,F8.6,8X,F8.6) 


GOTO 99 
ЕИ ҙ осо. 
C POLICY II (AGE REPLACEMENT - WITH MINIMAL REPAIR AND 
C SIMPLE PREVENTIVE MAINTENANCE) 
к D D a S 


250 CONTINUE 


PRINT*, ENTER THE EXPECTED COST OF A SIMPLE PREVENTIVE 
& MAINTENANCE FOLLOWED BY THE IMPROVEMENT IN THE 
& AGE IN HOURS ACHIEVED ON COMPLETION OF THE 
& MAINTENANCE’ 


READ*,C1P,XX 

PRINT*, EXPECTED COST OF SIMPLE PREVENTIVE MAINTENANCE 
& = ’,C1P, 

& "IMPROVEMENT IN AGE(HOURS) = „ХХ 


PRINT*, ENTER THE NUMBER OF SIMPLE PREVENTIVE 
& MAINTENANCE ACTIONS TO BE TAKEN IN A CYCLE ' 


READ*,NN 
PRINT*, NUMBER OF SIMPLE PREVENTIVE MAINTENANCE ACTIONS 
& = ММ 


PRINT*, ENTER THE VALUE FOR LOWER LIMIT, UPPER LIMIT, 
& STEP SIZE’ 

READ*,LL,UL,STEP 

PRINT*, LOWER LIMIT = ’,LL, UPPER LIMIT =’,UL, STEP = 

& "STEP 


WRITE(30,107)ALPHA,MU,C2P,C1C,C1P,XX,NN 
107 FORMAT(9X,POLICY,2X,1P,2X, (AGE REPLACEMENT WITH 
€ MINIMAL REPAIR AND SIMPLE PM) 
& ’/,9X,65(-’),//, 
& 10X,,ALPHA = ’,F3.1,17X,MU = ’,F6.1/, 
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& 10X,'C2P $,F7.1,12X/C1C = $,F7.1/, 

& 10X,C1P = $,F7.1,12X,XX = ’F6.1, HOURS’/, 

& 10X,"NO. OF SIMPLE PM =’,NN,//, 

& 5X, REPLACEMENT’,,6X,’COST’/, 5X, AGE (HOURS),5X,' 
& $/HOUR’ /,5%,110-),5%,7C-)) 


DO 6 J = LL,UL,STEP 
B = REAL(J) 
X = ХХ 
IF ((X-B) .GE. 0) THEN 
X=B 
END IF 


DO 7 I = 0,NN-1 
SUM1 = C1C*(LAMBDA**ALPHA)*(((B+(B-X)*D)**ALPHA) 
& - (((B-X)**ALPHA)*(I**ALPHA))) 
SUM = SUM + SUM1 
7 CONTINUE 


COST = (SUM + ((NN-1)*C1P) +C2P) / ((NN*B) +R2P) 
LENGTH = B*REAL(NN) 
WRITE(30,108) LENGTH,COST 
108 — FORMAT(SX,F6.1,5X,F9.4) 
SUM = 0.0 
6 CONTINUE 


GO TO 99 


C POLICY III (AGE REPLACEMENT - WITH TWO TYPES OF FAILURE) 


300 CONTINUE 


PRINT*, ENTER THE ESTIMATED VALUES FOR THE SHAPE 
& PARAMETER ALPHA AND THE SCALE PARAMETER MU’ 


READ*,ALPHA,MU 
PRINT*, ALPHA =’,ALPHA,’MU = ’,MU 


PRINT*, ENTER ESTIMATED VALUES FOR THE EXPECTED COST OF 
& PREVENTIVE REPLACEMENT FOLLOWED BY THE EXPECTED 
& COST OF FAILURE REPLACEMENT, AND COST OF MINIMAL 


& REPAIR’ 
READ*,C2P,C2C,C1C 


PRINT*, EXPECTED COST OF A PREVENTIVE REPLACEMENT = 
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',‚C2P 

‚EXPECTED COST OF A FAILURE REPLACEMENT = 
Доре: 

,» EXPECTED COST OF A MINIMAL REPAIR = 
СІС 


RINT*,ENTER VALUES FOR THE EXPECTED DOWNTIME OF 
PREVENTIVE REPLACEMENT FOLLOWED BY EXPECTED 
DOWNTIME OF A FAILURE REPLACEMENT AND DOWNTIME 
OF MINIMAL REPAIR’ 
READ*,R2P,R2C,R1C 
PRINT*'EXPECTED DOWNTIME OF A PREVENTIVE REPLACEMENT = 
',R2P 

, EXPECTED DOWNTIME OF A FAILURE REPLACEMENT = 
В2С 

, EXPECTED DOWNTIME OF A MINIMAL REPAIR = 
КІС 


€ gf uy 6060060 


Qo Ro Ro Ro Ro 


PRINT* ENTER THE VALUE FOR LOWER LIMIT, UPPER LIMIT, 
& STEP SIZE ' 

READ*,LL,UL,STEP 

PRINT*, LOWER LIMIT = ’,LL, UPPER LIMIT =’,UL,STEP = 

& "STEP 


PRINT*,/ENTER VALUES FOR THE PROBABILITY OF MINIMAL 
& FAILURE AND MAJOR FAILURE’ 

READ*,P1,P2 

PRINT*, THE PROBABILITY OF MINIMAL FAILURE = РЁ 
& AND PROBABILITY OF MAJOR FAILURE - x 


PRINT*/ENTER THE VALUE FOR LOWER LIMIT, UPPER LIMIT, 
& STEP SIZE ’ 

READ*,LL,UL,STEP 

PRINT*,, LOWER LIMIT = ’,LL, UPPER LIMIT =’,UL, STEP = 

& "STEP 


LAMBDA = 1.0D0/MU 


WRITE(30,109)ALPHA,MU,C2P,C2C,C1C,R2P,R2C,R1C,P1,P2 
109 FORMAT(9X,' POLICY',2X,'IIT',2X, (AGE REPLACEMENT WITH 
TWO TYPES OF FAILURE/,9X,55(-)//, 
$ 5X/ALPHA = ’,F3.1,11X,MU = ’,F6.1/, 
5X,C2P = $,F7.1,6X,C2C = $,F7.1,7X,C1C = 
$,F7.1/, SX,R2P = ',F4.1, HOURS’, 5X, R2C = 
',F4.1, HOURS’,5X, "RIC = ’,F4.1, HOURS’/,5X, 
"PROBABILITY OF TYPE I FAILURE = ’,F3.1,/,5X, 


Qo Ro Ro Ro Ro Ro 
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& PROBABILITY OF TYPE II FAILURE = ’,F3.1,//, 

& ,5X,REPLACEMENT’,6X, COST’,10X, APPROXIMATH’,8X, 
& ’EXACT’/, 5X, AGE (HOURS) ,5X,’$/HOUR’, 9X,’ 

& AVAILABILITY ’,4X, AVAILABILITY’ /,5X,11(-’), 

& 4Х,8(-),8Х,12(-), € 4X,12(-)) 


DO 13 Iz 0,4000,50 
B=DBLE(I) 


C CALL IMSL SUBROUTINE DQDAG 
CALL DQDAG (F,A,B,ERRABS,ERRREL,IRULE,RESULT,ERREST) 


А1 = DEXP(-((((LAMBDA*B)**ALPHA)*P1)+(B/R1C))) 
AA = Al+(RESULT*(1.0D0/R1C) ) 


C OUTPUT THE AVAILABILITY AT TIME T TO FILE 13 


WRITE(13,29)B,AA 
29 FORMAT(1X,F6.1,15X,F8.6) 


13 CONTINUE 
REWIND (13) 


DO 43 J = 1, NDATA 
READ(13,*)XDATA(J), YDATA(J) 
43 CONTINUE 


CALL DCSINT(NDATA,XDATA, YDATA,BREAK,CSCOEF) 


Calculate the integral of the spline 
approximation. 
DO 33 K = LL,UL,STEP 
B = DBLE(K) 
NINTV = NDATA - 1 
CALL DQDAG (G,A,B,ERRABS,ERRREL,IRULE,SOLN,ERREST) 
SF = DEXP(-(P2*((LAMBDA*B)**ALPHA))) 
CDF = 1.0D0-SF 
LENGTH = (R2P*SF)+(R2C*CDF)+SOLN 
AVAIL1 = SOLN/LENGTH 


DISS) 


ANSWER = DCSITG(A,B,NINTV,BREAK,CSCOEF) 
AVAIL = ANSWER/LENGTH 


CALL DQDAG (H,A,B,ERRABS,ERRREL,IRULE,ANS,ERREST) 
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COST = (C1C*SF*(P1*(LAMBDA*B)**ALPHA)) + (C2P*SF) + 
& (C2C*CDF) + (ALPHA*P1*P2*C1C*(LAMBDA** 
& (2.0DO*ALPHA)) *ANS) 


COSMIN = COST / LENGTH 


WRITE (30,110) B,COSMIN,AVAIL1,AVAIL 
33 CONTINUE 
110 FORMAT(7X,F6.1,6X,F9.5,10X,F8.6,8X,F8.6) 


99 STOP 
END 


DOUBLE PRECISION FUNCTION K(T) 
REAL*8 LAMBDA,ALPHA,T 

COMMON /HT/ LAMBDA,ALPHA 
E=DEXP(-((LAMBDA * T)**ALPHA)) 
RETURN 

END 


DOUBLE PRECISION FUNCTION EE(T) 

REAL*8 LAMBDA,ALPHA,T,R1C,B 

COMMON /HT/ LAMBDA,ALPHA 

COMMON /FT/ R1C,B 

EEzDEXP(-X(( LAMBDA * B)*FALPHA)-(LAMBDA*T)**ALPHA) + 
& (B/R1C) 

&-(T/R1C))) 

RETURN 

END 


DOUBLE PRECISION FUNCTION F(T) 
REAL*8 LAMBDA,ALPHA,T,R1C,B,P1,P2 
COMMON/FT/ R1C,B 

COMMON /HT/ LAMBDA,ALPHA 
COMMONV/IT/ P1 

COMMON /GT/ P2 
F=DEXP((((LAMBDA * B)**ALPHA)*P1) - 
& ((LAMBDA*T)**ALPHA)*P1)+(B/R1C) -(T/R1C) + 
& ((LAMBDA*T)**ALPHA)*P2) )) 
RETURN 

END 
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DOUBLE PRECISION FUNCTION FF(T) 

REAL*8 LAMBDA,ALPHA,T 

COMMON /HT/ LAMBDA,ALPHA 
FF=(T**ALPHA)*(DEXP(-((LAMBDA * T)**ALPHA))) 
RETURN 

END 


DOUBLE PRECISION FUNCTION G(T) 
REAL*8 LAMBDA,ALPHA,T,P2 
COMMON/GT/ P2 

COMMON /HT/ LAMBDA,ALPHA 
G=DEXP(-(P2*((LAMBDA*T)**ALPHA))) 
RETURN 

END 


DOUBLE PRECISION FUNCTION H(T) 

REAL*8 LAMBDA,ALPHA,T,P2 

COMMON/GT/ P2 

COMMON /HT/ LAMBDA,ALPHA 
=(T**((2.0D0*ALPHA)-1.0D0)) * 

& DEXP(-(P2*((LAMBDA*T)**ALPHA))) 

RETURN 

END 
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POLICY 1 (AGE REPLACEMENT) 


ALPHA = 3.0 MU = 1390.0 
C2P = $25000.0 C2C = $37500.0 
R2P = 8.0 HOURS R2C = 16.0 HOURS 


MISSION DURATION = 24.0 HOURS 


REPLACEMENT COST RATE AVAILABILITY MISSION RELIABILITY 
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AGE (HOURS) ($/HOUR) (A) (R) 
900.0 32.781 0.988396 0.977947 
950.0 31.900 0.988565 0.975493 

1000.0 31.178 0.988681 0.972916 
1050.0 30.592 0.988752 0.970214 
1100.0 30.123 0.988785 0.967391 
1150.0 29.754 0.988786 0.964447 
1200.0 29.472 0.988761 0.961383 
1250.0 29.263 0.988713 0.958200 
1300.0 29.115 0.988649 0.954899 
1350.0 29.020 0.988571 0.951482 
1400.0 28.968 0.988483 0.947951 
1450.0 28.951 0.988389 0.944306 
1500.0 28.962 0.988290 0.940549 
1550.0 28.995 0.988191 0.936681 
1600.0 29.043 0.988093 0.932704 


POLICY 1 (AGE REPLACEMENT WITH SIMPLE PREVENTIVE MAINTENANCE) 


ALPHA 
R2P = 
RIP = 


NO. OF INTERVAL 


€ O0 -1 O» Oi» GO N = 


= 3.0 
= 8.0 HOURS 


1.0 HOUR 


REPLACEMENT AGE (HOURS) AVAILABILITY 


917.6 
1743.4 
2486.7 
3155.6 
3757.7 
4299.5 
4787.1 
5226.0 
0621.0 
0976.5 
6296.5 
6584.4 
6843.6 
7076.8 
7286.7 
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MU 
R2C 
R 


1390.0 


48.0 HOURS 


0.1 


0.979421 
0.982197 
0.983305 
0.983591 
0.983585 
0.983480 
0.983344 
0.983206 
0.983075 
0.982955 
0.982848 
0.982751 
0.982664 
0.982587 
0.982518 


POLICY II (AGE REPLACEMENT WITH MINIMAL REPAIR) 
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ALPHA = 3.0 MU = 1390.0 
C2P = $25000.0 C1C = $37500.0 
R2P = 8.0 HOURS R1C = 1.0 HOURS 
REPLACEMENT COST APPROXIMATE EXACT 
AGE (HOURS) /HOUR AVAILABILITY AVAILABILITY 
1200.0 40.6693 0.992845 0.992847 
1300.0 42.5667 0.993258 0.993261 
1400.0 44.9682 0.993593 0.993595 
1500.0 47.8289 0.993862 0.993865 
1600.0 51.1154 0.994076 0.994080 
1700.0 04.8018 0.994245 0.994249 
1800.0 58.8682 0.994374 0.994379 
1900.0 63.2987 0.994469 0.994474 
2000.0 68.0807 0.994532 0.994539 
2100.0 73.2039 0.994569 0.994576 
2200.0 78.6598 0.994581 0.994589 
2300.0 84.4415 0.994571 0.994580 
2400.0 90.5432 0.994540 0.994551 
2500.0 96.9601 0.994490 0.994503 
2600.0 103.6880 0.994423 0.994437 


POLICY III (AGE REPLACEMENT WITH TWO TYPES OF FAILURE 
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ALPHA - 3.0 MU = 1390.0 

C2P = $25000.0 C2C = $37500.0 C1C = $ 1000.0 

R2P = 8.0 HOURS R2C = 24.0 HOURS RIC = 8.0 HOURS 

PROBABILITY OF TYPE I FAILURE = 0.6 
PROBABILITY OF TYPE II FAILURE = 0.4 

REPLACEMENT COST APPROXIMATE EXACT 

AGE (HOURS) $/HOUR AVAILABILITY AVAILABILITY 
1018.0 27.30690 0.989570 0.988997 
1042.0 26.88558 0.989629 0.989030 
1066.0 26.49157 0.989678 0.989053 
1090.0 26.12318 0.989718 0.989066 
1114.0 25.77884 0.989749 0.989070 
1138.0 25.45711 0.989772 0.989065 
1162.0 25.15666 0.989787 0.989052 
1186.0 24.87625 0.989794 0.989032 
1210.0 24.61474 0.989795 0.989004 
1234.0 24.37107 0.989789 0.988970 
1258.0 24.14424 0.989778 0.988929 
1282.0 23.93333 0.989760 0.988882 
1306.0 23.73747 0.989737 0.988830 
1330.0 23.55583 0.989709 0.988773 
1354.0 23.38767 0.989677 0.988710 
1378.0 23.23225 0.989639 0.988643 


APPENDIX B. DERIVATION OF COST RATE 


AND AVAILABILITY FUNCTIONS 
Let T be the time to failure 


СЕ 


T > t, 
Tue 
Then E[C] = C, P[T>t,] + C, P[TSt,] 


C, P[1 - F(t,)] * C; P[F(t,)] 


nn + Ta» t. 


{ R + T Е. 


Then Е [1] 


ЕП | T > t,] P[T > t,] + E[L | T < t,] P[T € t] 
Now, 


РО ТЕ o р ст 





S t] 
P[T S t,] 
( F (t) 
{ A 
= | F (t,) 
{ 1 ЕЕ 
Then, 


ar «БШ 


Therefore 


E[L] = 


(R, + t,) [1 - F(t,)] + Re (F(t) ЈРЕЛО 
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ERO СЕЕ) + R. [Fit,)] + t, [1 - F(t,)] + 


Ja t f(t) dt 


SEDI ЛЕ (ЕР ВК Е") *t,-t,F(t)l * 


ІК t fit) dt 


The integral |, ? t f(t) dt can be simplified by using 


| u dv = uv - Í v du 


t t 
then J t AD dt = +, F(t,) - рр ғо 
now E[L] = R, [1 - F(tp)] + R, [F(tp)] + tp- РЕФ 4 
= R, [1 — F(t,)] + R; [F(t,)] + Ја - № F(t) dt 


= в [1 - F(t] +R [F(t] + [* [1 - FQ) 4 


= R, [1 - F(t,)] + R, [F(t,)] + M(t,) 
where M(t,) = f. ? [1 - FO) dt 


E[C] 
Now C(t,) = 





E[L] 


с, РІ1- F(t,)] + C, P[F(t,)] 


B (1 - F(t,)] + Re [Ele] + J” - FO] d 
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Setting dC(t,) 





dt, 


we get C, £(t,) |н [| - F()) d - с #(+,) i [1 - F(0)) dt 


C, [F(t,)] - ©, [F(t,)] * C, IF(t))]^ - с. [F(t,)]° 
= C, - C, [F(t))] * £(t,) [R.C, - RC] 


rearranging the terms we get 


£t)  f? [1 - FO] dt С, - СЬ) - ЕЕ ( - F(t)) (, - €) 


C, {1 - F(t,)} + £(t,) [RC, — RC, 
= (=>) 


since h(t,) = 
1 - F(t,) 


then h(t,) m [| - F()) d (с, - с.) - Е(+,) (С, - с.) + 


с, * h(t) во, - ВС 
and this simplifies to 





C Ес, -К,С 
(t) f "t -FO)dr-FG = сте + Ве.) в 
Р P 


Now 
Mean life during a cycle 


Availability A(t,)= 
Expected length of cycle 
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t, [1 - FG] + J> t fit) dt 
R, (1 - FEDI + R, [FE] + t D - FEJ) + 7: Ло dt 
_ f. * [1 - F(t)] de 
R,[l - Ft) + R, [FG] + f. 20 - ғал 


Setting dA (t,) 





= 0 
dt 


ме вес fOO -FO dt + RLF (t,)] + Rll - F(t,)] } {1 - F(t,)} 


= j*p-F)d кіп-іІғеыі жығар - 6026001 1 


simplifying the equation above we get, 


ү” [1 - 2@)] dt + R [F (t.) ] + R, - R, [F (t,) ] E 
F (t,) N [1 - F@) dt - RIF(t)]? - R[F(t)] * RIEGO]? » 
t t t 
früu-rFola -rep f^n-ro& +1 fP0-FOja ) 


RIED] - 4 [f*[ -F0)d ) BI8G21 ) 


Cancelling and rearranging the terms we get 
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(ек) I E [1 - F] d 1 [R - R]} - { F(t,) [R; -R, - R[F (tp) ] 


+ R [F(.)] } = Rell - F(t,)] 


and this can be simplified to 


(ЕЕ I ЈР - FQ) dt ] (& - K1) - C F(&) IR -R]I2 -tF(t)] - 





R.[1 - F(tj)] 
Then 
Jü) јон. _ LN 
г до J” ü - F@ dt - FG) кк 
R, 





na) [7 U - FO de - FG) = ко 


Rearranging the terms we have 


t, K 1 
i [1 =- FO dt = (FO ED | 


R, - R, 
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APPENDIX C. DESCRIPTION OF IMSL SUBROUTINE DQDAG 


QDAG/DQDAG . (Single/Double precision) 


Purpose: 


Integrate a function using a globally adaptive scheme based on 


Gauss-Kronrod rules. 


Usage: 


Arguments 


Е - 


А Е 


5 Е 


CALL QDAG (F, A, B, ERRABS, ERRREL, IRULE, RESULT, 


ERREST) 


User-supplied FUNCTION to be integrated. The form is F(X), where 
X - Independent variable. (Input) 

F - The function value. (Output) 

F must be declared EXTERNAL in the calling program. 

Lower limit of integration. (Input) 


Upper limit of integration. (Input) 


ERRABS - Absolute accuracy desired. (Input) 


ERRREL - Relative accuracy desired. (Input) 


IRULE - 


Choice of quadrature rule. (Input) 
A Gauss-Kronrod rule is used with 
7 - 15 points if IRULE = 1 

10 - 21 points if IRULE = 2 


15 - 31 points if IRULE = 3 
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20 - 41 points if IRULE = 4 

25 - 51 points if IRULE = 5 

30 - 61 points if IRULE = 6 

IRULE = 2 is recommended for most functions. 

If the function has a peak singularity use IRULE = 1 

lf the function is oscillatory use IRULE = 6 
RESULT - Estimate of the integral from A to B of F. (Output) 


ERREST - Estimate of the absolute value of the error. (Output) 


Notes 

QDAG is a general-purpose integrator that uses a globally adaptive scheme 
in Order to reduce the absolute error. It subdivides the interval [A,B] and uses a 
(2k + 1)-point Gauss-Kronrod rule to estimate the integral over each subinterval. 
The error for each subinterval is estimated by comparison with the k-point Gauss 
quadrature rule. The subinterval with the largest estimated error is then bisected 
and the same procedure is applied to both halves. The bisection process is 
continued until either the error criterion is satisfied, roundoff error is detected, the 
subintervals become too small, or the maximum number of subintervals allowed 


is reached. 
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APPENDIX D. DESCRIPTION OF IMSL SUBROUTINE DCSINT 


CSINT/DCSINT 
Purpose: 


Usage: 


Arguments 


NDATA - 


XDATA - 


YDATA - 


BREAK - 


CSCOEF - 


Notes 


(Single/Double precision) 
Compute the cubic spline interpolant . 


CALL DCSINT (NDATA, XDATA, YDATA, BREAK, CSCOEF) 


Number of data points. (input) 

NDATA must be at least 2. 

Array of length NDATA containing the data point abscissas. 
(Input) 

Array of length NDATA containing the data point ordinates. 
(Input) 

Array of length NDATA containing the breakpoints for the 
piecewise cubic representation. (Output) 

Matnx of size 4 by NDATA containing the local coefficients of 


the cubic pieces. (Output) 


DCSINT computes the second denvative cubic spline interpolant to a set of 


data points (x,y,) for i = 1,2,...NDATA = N. The breakpoints of the spline are the 
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abscissas. Endpoint conditions are automatically determined by the program. 
These conditions correspond to the "not-a-knot" condition, which requires that the 


third derivative of the spline be continuous at the second and next-to-last 


breakpoint. 
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APPENDIX E. DESCRIPTION OF IMSL FUNCTION DCSITG 


CSITG/DCSITG  (Single/Double precision) 


Purpose: Evaluate the integral of a cubic spline 
Usage: CSITG(A, B, NINTV, BREAK, CSCOEF) 
Argument 
A - Lower limit of integration. (Input) 
B - Upper limit of integration. (Output) 
NINTV - Number of polynomial pieces. (Input) 
BREAK - Array of length NINTV+1 containing the breakpoints for the 


piecewise cubic representation. (Input) 
CSCOEF- Matrix of size 4 by NINTV+1 containing the local 
coefficients of the cubic pieces. (Input) 


DCSITG - Value of the integral of the spline from A to B. (Output) 


Notes 
DCSITG evaluates the integral of a cubic spline over an interval. A cubic 


spline is a piecewise polynomial of order 4. 
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